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Abstract 


The observable universe could be a 1+3-surface (the “brane”) embedded in a 1+3-+d- 
dimensional spacetime (the “bulk”), with Standard Model particles and fields trapped on the 
brane while gravity is free to access the bulk. At least one of the d extra spatial dimensions 
could be very large relative to the Planck scale, which lowers the fundamental gravity scale, 
possibly even down to the electroweak (~ TeV) level. This revolutionary picture arises in the 
framework of recent developments in M theory. The 1+10-dimensional M theory encompasses 
the known 1+9-dimensional superstring theories, and is widely considered to be a promising 
potential route to quantum gravity. At low energies, gravity is localized at the brane and 
general relativity is recovered, but at high energies gravity “leaks” into the bulk, behaving in 
a truly higher-dimensional way. This introduces significant changes to gravitational dynamics 
and perturbations, with interesting and potentially testable implications for high-energy astro- 
physics, black holes, and cosmology. Brane-world models offer a phenomenological way to test 
some of the novel predictions and corrections to general relativity that are implied by M theory. 
This review analyzes the geometry, dynamics and perturbations of simple brane-world mod- 
els for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane-worlds 
based on the Randall-Sundrum models. We also cover the simplest brane-world models in 
which 4-dimensional gravity on the brane is modified at low energies — the 5-dimensional 
Dvali—Gabadadze—Porrati models. Then we discuss co-dimension two branes in 6-dimensional 
models. 
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Update (21 April 2010) 

Inserted subsections on numerical solutions in Sections [6] and [7] and five figures. Added new 
Subsection [2.2] on the Randall-Sundrum model in string theory. Added new Section [9]on Dvali-— 
Gabadadze-—Porrati models and nine figures. Added new Section [10]on 6-dimensional models and 
one figure. More than 80 new references. 


Page 3: Extended a paragraph moved here from the end of Section 
Page 11: Added Subsection 2.2] “RS model in string theory”. 


Page 39: Added two paragraphs discussing the neglection of backreaction due to metric per- 
turbations in the bulk. 


Page 56: Added Subsection “Full numerical solutions”. 


Page 66: Removed former Figure 11: Damping of brane-world gravity waves on horizon re-entry 
due to massive mode generation. 


Page 67: Added Subsection “Full numerical solutions”. 
Page 76: Added Section [9] “DGP Models”. 
Page 87: Added Section [0] “6-Dimensional Models”. 


1 Introduction 


At high enough energies, Einstein’s theory of general relativity breaks down, and will be superceded 
by a quantum gravity theory. The classical singularities predicted by general relativity in gravi- 
tational collapse and in the hot big bang will be removed by quantum gravity. But even below 
the fundamental energy scale that marks the transition to quantum gravity, significant corrections 
to general relativity will arise. These corrections could have a major impact on the behaviour 
of gravitational collapse, black holes, and the early universe, and they could leave a trace — a 
“smoking gun” — in various observations and experiments. Thus it is important to estimate these 
corrections and develop tests for detecting them or ruling them out. In this way, quantum gravity 
can begin to be subject to testing by astrophysical and cosmological observations. 

Developing a quantum theory of gravity and a unified theory of all the forces and particles of 
nature are the two main goals of current work in fundamental physics. There is as yet no generally 
accepted (pre-)quantum gravity theory. Two of the main contenders are M theory (for reviews 
see, e.g., [359] [380]) and quantum geometry (loop quantum gravity; for reviews see, e.g., 
[412]). It is important to explore the astrophysical and cosmological predictions of both these 
approaches. This review considers only models that arise within the framework of M theory. 

In this review, we focus on RS brane-worlds (mainly the RS 1-brane model) and their gen- 
eralizations, with the emphasis on geometry and gravitational dynamics (see 
[351] 270, [363] [121] [49] [269] 272] for previous reviews with a broadly similar approach). Other 
reviews focus on string-theory aspects, e.g., [360], or on particle physics aspects, e.g., 
[75]. We also discuss the 5D DGP models, which modify general relativity 
at low energies, unlike the RS models; these models have become important examples in cos- 
mology for achieving late-time acceleration of the universe without dark energy. Finally, we give 
brief overviews of 6D models, in which the brane has co-dimension two, introducing very different 
features to the 5D case with co-dimension one branes. 


1.1 Heuristics of higher-dimensional gravity 


One of the fundamental aspects of string theory is the need for extra spatial dimensiond]. This 
revives the original higher-dimensional ideas of Kaluza and Klein in the 1920s, but in a new context 
of quantum gravity. An important consequence of extra dimensions is that the 

4-dimensional Planck scale Mp = M4 is no longer the fundamental scale, which is M44a, where 
d is the number of extra dimensions. This can be seen from the modification of the gravitational 
potential. For an Einstein—Hilbert gravitational action we have 


1 
Seravity = 5 pes dty —(4+d)g [aR = 244a : (1) 
2hay¢ 
1 
OG p= “ORs, - 5 “tdp CD ap = —A4+d Vo aB + Kad UHAT, p, (2) 
where X^ = (z#,yt,..., y), and Kj, q is the gravitational coupling constant, 
87 
Kaad = 81Gara = rad’ (3) 
Mīla 


The static weak field limit of the field equations leads to the 4 + d-dimensional Poisson equation, 
whose solution is the gravitational potential, 
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k 
V(r) x “itd : (4) 


1We do not consider timelike extra dimensions: see for an interesting example. 


If the length scale of the extra dimensions is L, then on scales r < L, the potential is 4 + d- 
dimensional, V ~ r~“+®, By contrast, on scales large relative to L, where the extra dimensions 
do not contribute to variations in the potential, V behaves like a 4-dimensional potential, i.e., r ~ L 
in the d extra dimensions, and V ~ L~¢r~!. This means that the usual Planck scale becomes an 
effective coupling constant, describing gravity on scales much larger than the extra dimensions, 
and related to the fundamental scale via the volume of the extra dimensions: 
MÈ ~ Mito LÀ. (5) 

If the extra-dimensional volume is Planck scale, i.e., L ~ My 1 then M4+a ~ Mp. But if the extra- 
dimensional volume is significantly above Planck scale, then the true fundamental scale M4+q can 
be much less than the effective scale M, ~ 10!° GeV. In this case, we understand the weakness of 
gravity as due to the fact that it “spreads” into extra dimensions and only a part of it is felt in 4 
dimensions. 

A lower limit on M444 is given by null results in table-top experiments to test for deviations from 
Newton’s law in 4 dimensions, V x r~!. These experiments currently probe sub-millimetre 
scales, so that 


L107! mm ~ (1071 TeV => Maşa Z 108 TeV. (6) 


Stronger bounds for brane-worlds with compact flat extra dimensions can be derived from null 
results in particle accelerators and in high-energy astrophysics [195]. 


1.2 Brane-worlds and M theory 


String theory thus incorporates the possibility that the fundamental scale is much less than the 
Planck scale felt in 4 dimensions. There are five distinct 1+9-dimensional superstring theories, 
all giving quantum theories of gravity. Discoveries in the mid-1990s of duality transformations 
that relate these superstring theories and the 1+10-dimensional supergravity theory, led to the 
conjecture that all of these theories arise as different limits of a single theory, which has come to 
be known as M theory. The 11th dimension in M theory is related to the string coupling strength; 
the size of this dimension grows as the coupling becomes strong. At low energies, M theory can be 
approximated by 1+10-dimensional supergravity. 

It was also discovered that p-branes, which are extended objects of higher dimension than strings 
(1-branes), play a fundamental role in the theory. In the weak coupling limit, p-branes (p > 1) 
become infinitely heavy, so that they do not appear in the perturbative theory. Of particular 
importance among p-branes are the D-branes, on which open strings can end. Roughly speaking, 
open strings, which describe the non-gravitational sector, are attached at their endpoints to branes, 
while the closed strings of the gravitational sector can move freely in the bulk. Classically, this is 
realised via the localization of matter and radiation fields on the brane, with gravity propagating 
in the bulk (see Figure [I}. 

In the Horava-Witten solution [204], gauge fields of the standard model are confined on two 
1+9-branes located at the end points of an $1/Z_ orbifold, i.e., a circle folded on itself across a 
diameter. The 6 extra dimensions on the branes are compactified on a very small scale close to 
the fundamental scale, and their effect on the dynamics is felt through “moduli” fields, i.e., 5D 
scalar fields. A 5D realization of the Horava—Witten theory and the corresponding brane-world 
cosmology is given in [802] [303] [304]. 

These solutions can be thought of as effectively 5-dimensional, with an extra dimension that can 
be large relative to the fundamental scale. They provide the basis for the Randall-Sundrum (RS) 
2-brane models of 5-dimensional gravity (see Figure). The single-brane Randall—Sundrum 
models [361] with infinite extra dimension arise when the orbifold radius tends to infinity. The RS 


Figure 1: Schematic of confinement of matter to the brane, while gravity propagates in the bulk 


(from [75]). 
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Figure 2: The RS 2-brane model. (Figure taken from .) 


models are not the only phenomenological realizations of M theory ideas. They were preceded by 
the Arkani-Hamed—Dimopoulos—Dvali (ADD) brane-world models [15] [14] [13] [6] [370) [424] (169) [774], 
which put forward the idea that a large volume for the compact extra dimensions would lower the 
fundamental Planck scale, 


Mew ~ 1 TeV S Maza < Mp ~ 1078 TeV, (7) 


where Mew is the electroweak scale. If M414 is close to the lower limit in Equation (7), then this 
would address the long-standing “hierarchy” problem, i.e., why there is such a large gap between 
Mew and Mp- 

In the ADD models, more than one extra dimension is required for agreement with experiments, 
and there is “democracy” amongst the equivalent extra dimensions, which are typically flat. By 
contrast, the RS models have a “preferred” extra dimension, with other extra dimensions treated 
as ignorable (i.e., stabilized except at energies near the fundamental scale). Furthermore, this extra 
dimension is curved or “warped” rather than flat: The bulk is a portion of anti-de Sitter (AdSs5) 
spacetime. As in the Horava—Witten solutions, the RS branes are Z2-symmetric (mirror symmetry), 
and have a tension, which serves to counter the influence of the negative bulk cosmological constant 
on the brane. This also means that the self-gravity of the branes is incorporated in the RS models. 
The novel feature of the RS models compared to previous higher-dimensional models is that the 
observable 3 dimensions are protected from the large extra dimension (at low energies) by curvature 
rather than straightforward compactification. 

The RS brane-worlds and their generalizations (to include matter on the brane, scalar fields 
in the bulk, etc.) provide phenomenological models that reflect at least some of the features of 
M theory, and that bring exciting new geometric and particle physics ideas into play. The RS 
models also provide a framework for exploring holographic ideas that have emerged in M theory. 
Roughly speaking, holography suggests that higher-dimensional gravitational dynamics may be 
determined from knowledge of the fields on a lower-dimensional boundary. The AdS/CFT corre- 
spondence is an example, in which the classical dynamics of the higher-dimensional gravitational 
field are equivalent to the quantum dynamics of a conformal field theory (CFT) on the boundary. 
The RS model with its AdS; metric satisfies this correspondence to lowest perturbative order [130] 
(see also for the AdS/CFT correspondence in a cosmological 
context). 

Before turning to a more detailed analysis of RS brane-worlds, We discuss the notion of Kaluza— 
Klein (KK) modes of the graviton. 


1.3 Heuristics of KK modes 


The dilution of gravity via extra dimensions not only weakens gravity on the brane, it also extends 
the range of graviton modes felt on the brane beyond the massless mode of 4-dimensional gravity. 
For simplicity, consider a flat brane with one flat extra dimension, compactified through the iden- 
tification y © y + 2mnL, where n = 0,1,2,.... The perturbative 5D graviton amplitude can be 
Fourier expanded as 

fe w=) Fe"); (8) 

n 

where f, are the amplitudes of the KK modes, i.e., the effective 4D modes of the 5D graviton. 
To see that these KK modes are massive from the brane viewpoint, we start from the 5D wave 
equation that the massless 5D field f satisfies (in a suitable gauge): 


Oof=0 => Of +a7f=o. (9) 
It follows that the KK modes satisfy a 4D Klein—Gordon equation with an effective 4D mass Mn, 


— m2 = 
Ín = Ma fn, Mn = L (10) 


The massless mode fo is the usual 4D graviton mode. But there is a tower of massive modes, 
L~',2L~',..., which imprint the effect of the 5D gravitational field on the 4D brane. Compactness 
of the extra dimension leads to discreteness of the spectrum. For an infinite extra dimension, 
L —> œ, the separation between the modes disappears and the tower forms a continuous spectrum. 
In this case, the coupling of the KK modes to matter must be very weak in order to avoid exciting 
the lightest massive modes with m = 0. 

From a geometric viewpoint, the KK modes can also be understood via the fact that the 
projection of the null graviton 5-momentum p4 onto the brane is timelike. If the unit normal to 
the brane is n4, then the induced metric on the brane is 


) 


JAB = g JAB — NANB, Ogapnan® =1, gaBn” =0, (11) 


and the 5-momentum may be decomposed as 
Op 4 = MNA + PA, pan^ = 0, m= Op 4 nA, (12) 


where pa = gap p® is the projection along the brane, depending on the orientation of the 
5-momentum relative to the brane. The effective 4-momentum of the 5D graviton is thus p4. 
Expanding gap p4 ©)p® = 0, we find that 


gaBp p? = —m’. (13) 


It follows that the 5D graviton has an effective mass m on the brane. The usual 4D graviton 
corresponds to the zero mode, m = 0, when “p, is tangent to the brane. 

The extra dimensions lead to new scalar and vector degrees of freedom on the brane. In 5D, 
the spin-2 graviton is represented by a metric perturbation h4, that is transverse traceless: 


(pA, =0= dp Ohg®. (14) 


In a suitable gauge, hAg contains a 3D transverse traceless perturbation hij, a 3D transverse 
vector perturbation 4;, and a scalar perturbation 6, which each satisfy the 5D wave equation [131]: 


Chan — hij, Da, 8, (15) 
hi; = 0 = hË, (16) 
8,5! =0, (17) 

B 
(O+ 67) | = | =0. (18) 


The other components of hap are determined via constraints once these wave equations are 
solved. The 5 degrees of freedom (polarizations) in the 5D graviton are thus split into 2 (hi;) + 2 
(5;) +1 (8) degrees of freedom in 4D. On the brane, the 5D graviton field is felt as 


e a 4D spin-2 graviton hij (2 polarizations), 
e a 4D spin-1 gravi-vector (gravi-photon) ©; (2 polarizations), and 
e a 4D spin-0 gravi-scalar 8. 


The massive modes of the 5D graviton are represented via massive modes in all three of these fields 
on the brane. The standard 4D graviton corresponds to the massless zero-mode of hij. 

In the general case of d extra dimensions, the number of degrees of freedom in the graviton 
follows from the irreducible tensor representations of the isometry group as 4(d + 1)(d+ 4). 


2 Randall—Sundrum Brane-Worlds 


RS brane-worlds do not rely on compactification to localize gravity at the brane, but on the 
curvature of the bulk (sometimes called “warped compactification” ). What prevents gravity from 
‘leaking’ into the extra dimension at low energies is a negative bulk cosmological constant, 


-z= (19) 
where £ is the curvature radius of AdS; and y is the corresponding energy scale. The curvature 
radius determines the magnitude of the Riemann tensor: 

(5) 


(5) (5) 


1 
Rascp = -7 [aac gap — gan gael . (20) 
The bulk cosmological constant acts to “squeeze” the gravitational field closer to the brane. We 
can see this clearly in Gaussian normal coordinates X4 = (2, y) based on the brane at y = 0, for 


which the AdS; metric takes the form 
(ds? = e?n dada” + dy’, (21) 


with jy, being the Minkowski metric. The exponential warp factor reflects the confining role of 
the bulk cosmological constant. The Zo-symmetry about the brane at y = 0 is incorporated via 
the |y| term. In the bulk, this metric is a solution of the 5D Einstein equations, 


OGap=—As gar, (22) 


i.e., © T4_ = 0in Equation (2). The brane is a flat Minkowski spacetime, g4p(x%,0) = Nuvd" Ao” B, 
with self-gravity in the form of brane tension. One can also use Poincare coordinates, which bring 
the metric into manifestly conformally flat form, 


(2 
ds? = < [nu dz" dz" + d2?°], (23) 
zZ 
where z = le¥/*, 


The two RS models are distinguished as follows: 


RS 2-brane: There are two branes in this model [362], at y = 0 and y = L, with Z-symmetry 
identifications 
yor-y, yt+tLoL—y. (24) 


The branes have equal and opposite tensions +A, where 


3M2 
= EE (25) 


The positive-tension brane has fundamental scale Ms and is “hidden”. Standard model fields 
are confined on the negative tension (or “visible” ) brane. Because of the exponential warping 
factor, the effective Planck scale on the visible brane at y = L is given by 


M? = MŠ? Jerre = 1| . (26) 


So the RS 2-brane model gives a new approach to the hierarchy problem: even if Ms ~ 47t ~ 
TeV, we can recover Mp ~ 1016 TeV by choosing L/£ large enough. Because of the finite 


separation between the branes, the KK spectrum is discrete. Furthermore, at low energies 
gravity on the branes becomes Brans—Dicke-like, with the sign of the Brans—Dicke parameter 
equal to the sign of the brane tension [156]. In order to recover 4D general relativity at low 
energies, a mechanism is required to stabilize the inter-brane distance, which corresponds to 


a scalar field degree of freedom known as the radion [285]. 


RS 1-brane: In this model [361], there is only one, positive tension, brane. It may be thought of 
as arising from sending the negative tension brane off to infinity, L —> oo. Then the energy 


scales are related via 
MÈ 


£ 


The infinite extra dimension makes a finite contribution to the 5D volume because of the 


warp factor: 
[EXV alas] dye™*”/t = 5 f ate. (28) 
0 


Thus the effective size of the extra dimension probed by the 5D graviton is £. 


M3 = 


(27) 


We will concentrate mainly on RS 1-brane from now on, referring to RS 2-brane occasionally. The 
RS 1-brane models are in some sense the most simple and geometrically appealing form of a brane- 
world model, while at the same time providing a framework for the AdS/CFT correspondence [130 
[81]. The RS 2-brane introduce the added complication of radion 
stabilization, as well as possible complications arising from negative tension. However, they remain 
important and will occasionally be discussed. 


2.1 KK modes in RS 1-brane 


In RS 1-brane, the negative As is offset by the positive brane tension A. The fine-tuning in 
Equation ensures that there is a zero effective cosmological constant on the brane, so that the 
brane has the induced geometry of Minkowski spacetime. To see how gravity is localized at low 
energies, we consider the 5D graviton perturbations of the metric 123], 


(5) ) 


gap + e7211 (5) 


JAB `> (5 hap, Oh ay =0= Dh = Spr p (29) 


(see Figure B). This is the RS gauge, which is different from the gauge used in Equation (15), but 
which also has no remaining gauge freedom. The 5 polarizations of the 5D graviton are contained 
in the 5 independent components of Ohy in the RS gauge. 

We split the amplitude f of ©h4pg into 3D Fourier modes, and the linearized 5D Einstein 
equations lead to the wave equation (y > 0) 


HKF =e r — ; r'| . (30) 


Separability means that we can write 


Fy) =X emit) fm lU), (31) 


and the wave equation reduces to 
Pm + (m? + Kk?) om = 0, (32) 
n S Sin +e? fn = 0. (33) 


Figure 3: Gravitational field of a small point particle on the brane in RS gauge. (Figure taken 
from [156].) 


The zero mode solution is 


polt) = Aope tt + Ape, (34) 
folu) = Bo + Coet”, (35) 
and the m > 0 solutions are 
Ym (t) = Am4 exp ( hiy m? +k? t) + Am- exp (-i m2 +k? t) , (36) 
fm (y) = e/t [Bm Jo (meet/*) + Co (meen!) ] (37) 


The boundary condition for the perturbations arises from the junction conditions, Equa- 
tion (68), discussed below, and leads to f’(t,0) = 0, since the transverse traceless part of the 
perturbed energy-momentum tensor on the brane vanishes. This implies 


_ Ii(me) 


= m= Br: 


The zero mode is normalizable, since 


J Boe?!" < 00. (39) 
0 


Its contribution to the gravitational potential V = 4 (hog gives the 4D result, V œ r~!. The 
contribution of the massive KK modes sums to a correction of the 4D potential. For r < £Z, one 


obtains 
GMe 


2 2 


(40) 


(ne 


which simply reflects the fact that the potential becomes truly 5D on small scales. For r > 4, 


V(r) & cu (1 + a) ; (41) 


3r2 


10 


which gives the small correction to 4D gravity at low energies from extra-dimensional effects. These 
effects serve to slightly strengthen the gravitational field, as expected. 

Table-top tests of Newton’s laws currently find no deviations down to O(10~' mm), so that 
L 0.1 mm in Equation (41). Then by Equations and (27), this leads to lower limits on the 
brane tension and the fundamental scale of the RS 1-brane model: 


à > (1 TeV)*, M; > 10° TeV. (42) 


These limits do not apply to the 2-brane case. 

For the 1-brane model, the boundary condition, Equation (88), admits a continuous spectrum 
m > 0 of KK modes. In the 2-brane model, f’(t, L) = 0 must hold in addition to Equation (38). 
This leads to conditions on m, so that the KK spectrum is discrete: 


Jy (me) 
Yı (me) 


In _ 
Tn ,-L/e 


Mn = where Ji (£n) = Yi (Tn): (43) 
The limit Equation (42) indicates that there are no observable collider, i.e., O(TeV), signatures 
for the RS 1-brane model. The 2-brane model by contrast, for suitable choice of L and £ so 
that mı = O(TeV), does predict collider signatures that are distinct from those of the ADD 


models 195). 


2.2 RS model in string theory 


There have been interesting developments in realizing the Randall-Sundrum brane-world model in 
string theory. A concrete example was found in type IIB supergravity [171]. The bosonic part of 
the 10-dimensional action is given by (see for example [358]) 


acy x 1 io VIJ —2¢ [R A(V 2) Fa ri G G Fé 
cee a Oi) A es a ae 
1 = 
toe | Cw A Go) A Gos) + Soc. (44) 
iko 


Here gs denotes the string metric. We have also defined the combined three-flux, G3) = F(3) — 
TH (3), where T = Co) + ie~?, and 

~ 1 1 

Fis) = Fis) — 5 Ca) ^ Hay + 5B ^ Fo). (45) 


The term Sloc is the action of localized objects, such as branes. 
There is a 10-dimensional solution similar to the one in the RS model 


ds?o = h(y)~? mvda"da” + h(y)? Gnndy™dy”, (46) 


where x" are four-dimensional coordinates and y™ are coordinates on the 6-dimensional compact 
manifold Mg. If there are N coincident D3-branes on the manifold, the 10-dimensional spacetime 
is described by AdSs x Xs where X; is a five-dimensional Einstein manifold. The warp factor is 


given by 
4 
rom (=) (47) 


where r is the distance from the D3 branes in the gmn metric and L is given by 


L =4nag,Na”. (48) 
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Here gs is the string coupling constant, a’ is related to the string mass scale ms as a’ = 1/m? and 
a depends on Xs. 

Near the D3 branes, the geometry is AdS x S5. At large r, the product structure of the 
10-dimensional metric breaks down and r ceases to be a good coordinate, and the 6D extra- 
dimensional space Mg is not a product of a five sphere and one-dimensional spacetime. This gives 
a minimum distance ro. This cut-off acts as a TeV brane in the RS model. On the other hand, 
for r > Tmax, the AdS geometry smoothly glues into a Calabi-Yau compactification. Thus, the 
compact manifold plays the role of the Planck brane in the RS model, acting as a cut-off of the 
spacetime. The minimum warp factor h(r9)~!/? is determined by the flux associated with F3 and 
H; and thanks to the warping, it is possible to realize a large hierarchy 


h(ro) 


irony 2h (49) 
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3  Covariant Approach to Brane-World Geometry and Dy- 
namics 


The RS models and the subsequent generalization from a Minkowski brane to a Friedmann- 
Robertson—Walker (FRW) brane were derived as solutions 
in particular coordinates of the 5D Einstein equations, together with the junction conditions at 
the Zo-symmetric brand}. A broader perspective, with useful insights into the inter-play between 
4D and 5D effects, can be obtained via the covariant Shiromizu-Maeda-Sasaki approach [391], in 
which the brane and bulk metrics remain general. The basic idea is to use the Gauss—Codazzi 
equations to project the 5D curvature along the brane. (The general formalism for relating the 
geometries of a spacetime and of hypersurfaces within that spacetime is given in [426].) 
The 5D field equations determine the 5D curvature tensor; in the bulk, they are 


OG ag =—As gag + K? Tap, (50) 


where )T4pg represents any 5D energy-momentum of the gravitational sector (e.g., dilaton and 
moduli scalar fields, form fields). 

Let y be a Gaussian normal coordinate orthogonal to the brane (which is at y = 0 without loss 
of generality), so that n4dX4 = dy, with n^ being the unit normal. The 5D metric in terms of 
the induced metric on {y = const.} surfaces is locally given by 


YaB = JAB t NANB, Ods? = Gould”; y)dax! da” + dy’. (51) 


The extrinsic curvature of {y = const. } surfaces describes the embedding of these surfaces. It can 
be defined via the Lie derivative or via the covariant derivative: 


1 
KaB = gin gap = ga° Von, (52) 


so that 
Kian) =0= Kapn”, (53) 


where square brackets denote anti-symmetrization. The Gauss equation gives the 4D curvature 
tensor in terms of the projection of the 5D curvature, with extrinsic curvature corrections: 


Rapop = Rercuga’ gB” 90° 9p" + 2K acK pp, (54) 
and the Codazzi equation determines the change of K4pg along {y = const.} via 
Vek? 4- VAK = Rao gab”, (55) 


where K = K4 4. 
Some other useful projections of the 5D curvature are: 


ORerau ga” gB” goln” = 2V4 Kpc, (56) 
ORprag gafn" ggn” = —-L£y)K apt Kack® sz, (57) 
ORop gac gp? = Rap — £nKap—KKapt2Kack® p. (58) 


The 5D curvature tensor has Weyl (tracefree) and Ricci parts: 


2 1 
ORapop = © Cacep + 3 (Osace ORD Bp — gzic Rp) = o gale Ogn OR. (59) 


2Brane-worlds without Z2 symmetry have also been considered; see e.g., for recent work. 
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3.1 Field equations on the brane 


Using Equations (50) and (64), it follows that 


Gw = — 5 Msg + ond OTa Bgu g” + (ETag? — L er) ow 
+KK pw K, Kov + ; [K Kap — K?| gw — Ew, (60) 
where ÖT = )T4 4, and where 
Ew = OCacep ae ae ae (61) 


is the projection of the bulk Weyl tensor orthogonal to nê. This tensor satisfies 
ann” =0 = Ejaan) =a", (62) 


by virtue of the Weyl tensor symmetries. Evaluating Equation (60) on the brane (strictly, as 
y — +0, since Eag is not defined on the brane [391]) will give the field equations on the brane. 

First, we need to determine K, at the brane from the junction conditions. The total energy- 
momentum tensor on the brane is 


INS = Tuv = AG (63) 


where Tuv is the energy-momentum tensor of particles and fields confined to the brane (so that 
Tapn? =0). The 5D field equations, including explicitly the contribution of the brane, are then 


OG ap =—As gap + Kk? [Tas + T38” 5(y)| . (64) 


Here the delta function enforces in the classical theory the string theory idea that Standard Model 
fields are confined to the brane. This is not a gravitational confinement, since there is in general 
a nonzero acceleration of particles normal to the brane [805]. 

Integrating Equation (64) along the extra dimension from y = —e to y = +e, and taking the 
limit € + 0, leads to the Israel—-Darmois junction conditions at the brane, 


Gav — Suv = 9, (65) 


none ees 
Ky Ko = Ake |e = P| (66) 


where Toe? = grr, The Zə symmetry means that when you approach the brane from one 
side and go through it, you emerge into a bulk that looks the same, but with the normal reversed, 
nê — —n4. Then Equation implies that 


Ki =—Kj}, (67) 


uv? 


so that we can use the junction condition Equation to determine the extrinsic curvature on 
the brane: 


1 1 
Kuw = -3545 Tu +3 (AT) Suv ; (68) 


where T = T“,,, where we have dropped the (+), and where we evaluate quantities on the brane 
by taking the limit y + +0. 
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Finally we arrive at the induced field equations on the brane, by substituting Equation 
into Equation (60): 
Pe. 2 


Gu = Agw H K Tiv H 6 Su T Cu Ea 4 Fy (69) 


The 4D gravitational constant is an effective coupling constant inherited from the fundamental 
coupling constant, and the 4D cosmological constant is nonzero when the RS balance between the 
bulk cosmological constant and the brane tension is broken: 

1 


Kk = res (70) 


A= 5 [As +A]. (71) 


The first correction term relative to Einstein’s theory is quadratic in the energy-momentum 
tensor, arising from the extrinsic curvature terms in the projected Einstein tensor: 
1 1 
TT yy — i 


Sw = 7p 


1 
TyoT*y + sy 9uv [3Tas T? — T°]. (72) 


The second correction term is the projected Weyl term. The last correction term on the right of 
Equation (69), which generalizes the field equations in [391], is 


l; 


F ww = OTa Bgu g” F OT sapn4n® A a Juv, (73) 


where Č)T4pg describes any stresses in the bulk apart from the cosmological constant (see [313] for 
the case of a scalar field). 
What about the conservation equations? Using Equations (50), and (68), one obtains 


WT = —2 igh 6". (74) 


Thus in general there is exchange of energy-momentum between the bulk and the brane. From 
now on, we will assume that 


Tag 0 > Fi =O, (75) 
so that 
OGap =—As gan in the bulk, (76) 
od 
Cin Ag Pe Ti 6S — Euv on the brane. (77) 


One then recovers from Equation (74) the standard 4D conservation equations, 
V Tush: (78) 


This means that there is no exchange of energy-momentum between the bulk and the brane; their 
interaction is purely gravitational. Then the 4D contracted Bianchi identities (V”G „v = 0), applied 
to Equation (69), lead to 

6K? 

V” Ew = ` V“ Suv, (79) 
which shows qualitatively how 1 + 3 spacetime variations in the matter-radiation on the brane can 
source KK modes. 

The induced field equations (77) show two key modifications to the standard 4D Einstein field 
equations arising from extra-dimensional effects: 
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© Sw ~ (Ty)? is the high-energy correction term, which is negligible for p « A, but dominant 
for p > à: 
|k?Su/A| > [Zw P (80) 
\K?Tv| À A 
e €,, is the projection of the bulk Weyl tensor on the brane, and encodes corrections from 5D 
graviton effects (the KK modes in the linearized case). From the brane-observer viewpoint, 
the energy-momentum corrections in S,, are local, whereas the KK corrections in Epy are 
nonlocal, since they incorporate 5D gravity wave modes. These nonlocal corrections cannot 
be determined purely from data on the brane. In the perturbative analysis of RS 1-brane 
which leads to the corrections in the gravitational potential, Equation (41), the KK modes 
that generate this correction are responsible for a nonzero €,,; this term is what carries the 
modification to the weak-field field equations. The 9 independent components in the tracefree 
Env are reduced to 5 degrees of freedom by Equation (79); these arise from the 5 polarizations 
of the 5D graviton. Note that the covariant formalism applies also to the two-brane case. In 
that case, the gravitational influence of the second brane is felt via its contribution to Ev. 


3.2 5-dimensional equations and the initial-value problem 


The effective field equations are not a closed system. One needs to supplement them by 5D 
equations governing €,,,, which are obtained from the 5D Einstein and Bianchi equations. This 
leads to the coupled system [877] 


1 
En Ky = KaR y = Eno T gi5Im» (81) 
Q 1 7 Q Q 
EnEn = V" Bana + gs (Kuv — QuvK) + K” Ruavp + 8K (wEvja — KEw 
+ (Kyakvp — KopK yw) K”, (82) 
LnBuve = —2V pEvja alr Kaf Buvp — 2B ala (83) 
Ln Ruvag = —2RyvyfaK gy" a VuBapy T VpBpav, (84) 


where the “magnetic” part of the bulk Weyl tensor, counterpart to the “electric” part Ey, is 
Buva = gu gv? ga? © Cascpn”. (85) 
These equations are to be solved subject to the boundary conditions at the brane, 
Voi = KEV" Spv, (86) 
Buva = 2V u Kija = -Viu (Te = sive (87) 
where A = B denotes Alprane = B|brane- 
The above equations have been used to develop a covariant analysis of the weak field [377]. 


They can also be used to develop a Taylor expansion of the metric about the brane. In Gaussian 
normal coordinates, Equation (61), we have £n = 0/Oy. Then we find 


al 
E EE =e z POs Tjauw | 
3 y=0+ 


1 2 1 
+ [es + z's Ge + ri a Ty) + 6 (Gato i T) ~ As) ow y? F Sse 
Yy 
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In a non-covariant approach based on a specific form of the bulk metric in particular coordinates, 
the 5D Bianchi identities would be avoided and the equivalent problem would be one of solving 
the 5D field equations, subject to suitable 5D initial conditions and to the boundary conditions 
Equation on the metric. The advantage of the covariant splitting of the field equations 
and Bianchi identities along and normal to the brane is the clear insight that it gives into the 
interplay between the 4D and 5D gravitational fields. The disadvantage is that the splitting is 
not well suited to dynamical evolution of the equations. Evolution off the timelike brane in the 
spacelike normal direction does not in general constitute a well-defined initial value problem [I2]. 
One needs to specify initial data on a 4D spacelike (or null) surface, with boundary conditions 
at the brane(s) ensuring a consistent evolution [208]. Clearly the evolution of the observed 
universe is dependent upon initial conditions which are inaccessible to brane-bound observers; this 
is simply another aspect of the fact that the brane dynamics is not determined by 4D but by 
5D equations. The initial conditions on a 4D surface could arise from models for creation of the 
5D universe [382], from dynamical attractor behaviour [337] or from suitable 
conditions (such as no incoming gravitational radiation) at the past Cauchy horizon if the bulk is 
asymptotically AdS. 


3.3 The brane viewpoint: A 1-+ 3-covariant analysis 


Following |305), a systematic analysis can be developed from the viewpoint of a brane-bound ob- 
server. (See also [229].) The effects of bulk gravity are conveyed, from a brane observer viewpoint, 
via the local (S,,,) and nonlocal (E y) corrections to Einstein’s equations. (In the more general 
case, bulk effects on the brane are also carried by F,,,, which describes any 5D fields.) The Epy 
term cannot in general be determined from data on the brane, and the 5D equations above (or 
their equivalent) need to be solved in order to find E,,,. 

The general form of the brane energy-momentum tensor for any matter fields (scalar fields, 
perfect fluids, kinetic gases, dissipative fluids, etc.), including a combination of different fields, can 
be covariantly given in terms of a chosen 4-velocity u” as 


Taw = Puppy + phy + Ty + pty + qup- (89) 
Here p and p are the energy density and isotropic pressure, respectively, and 


(5) 


Rue = Juv + yt = guv — Nate + Ugly (90) 


projects into the comoving rest space orthogonal to u” on the brane. The momentum density and 
anisotropic stress obey 


lu =Un)> Tuv = Myr), (91) 
where angled brackets denote the spatially projected, symmetric, and tracefree part: 


V Ex 1 a 
Viy =P Vor Win = [hhf — 5h? hal Wag. (92) 
In an inertial frame at any point on the brane, we have 
u” = (1,0), hyv = diag(0, 1,1,1), Va = (0, V:),  Wpo =0= X Wa = Wij—-Wji, (93) 


where 2,7 = 1,2,3. 
The tensor S,,, which carries local bulk effects onto the brane, may then be irreducibly de- 
composed as 


1 a 1 Q Q 
Si = om [2p” — 3napT | Upp + rl [2p + 4pp + Tagt Po 4qaq ] huv 
1 


p+ 3P)T uv — Talu Tv) + Fu) + 3PUuv) — 541 Ta(puy): (94) 


“Tp! 4 2 
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This simplifies for a perfect fluid or minimally-coupled scalar field to 


Suv = <p [Pupun + (P + 2p) hyr]: (95) 


l 
=P’ 

The tracefree €,, carries nonlocal bulk effects onto the brane, and contributes an effective 
“dark” radiative energy-momentum on the brane, with energy density pg, pressure pe/3, momen- 
tum density qË, and anisotropic stress mi 


1 1 
Ew = pe (inw + D) + Eu + 9° uy + ce (96) 
We can think of this as a KK or Weyl “fluid”. The brane “feels” the bulk gravitational field 
through this effective fluid. More specifically: 


e The KK (or Weyl) anisotropic stress m£, incorporates the scalar or spin-0 (“Coulomb”), the 


vector (transverse) or spin-1 (gravimagnetic), and the tensor (transverse traceless) or spin-2 
(gravitational wave) 4D modes of the spin-2 5D graviton. 


e The KK momentum density gi incorporates spin-O and spin-1 modes, and defines a velocity 
ve of the Weyl fluid relative to u via gi = pevs,. 


e The KK energy density pg, often called the “dark radiation”, incorporates the spin-0 mode. 


In special cases, symmetry will impose simplifications on this tensor. For example, it must 
vanish for a conformally flat bulk, including AdSs, 


(9g 4% conformally flat => Ew =0. (97) 


The RS models have a Minkowski brane in an AdSs bulk. This bulk is also compatible with an FRW 
brane. However, the most general vacuum bulk with a Friedmann brane is Schwarzschild-anti-de 
Sitter spacetime [45]. Then it follows from the FRW symmetries that 

Schwarzschild AdSsbulk, FRW brane: 4 = 0 = 1%, 


pv? 


(98) 


where pe = 0 only if the mass of the black hole in the bulk is zero. The presence of the bulk black 
hole generates via Coulomb effects the dark radiation on the brane. 
For a static spherically symmetric brane (e.g., the exterior of a static star or black hole) [105], 


static spherical brane: a. = 0. (99) 


This condition also holds for a Bianchi I brane [309]. In these cases, ites is not determined by 


the symmetries, but by the 5D field equations. By contrast, the symmetries of a Gödel brane fix 
E 
Taw [25]. 
The brane-world corrections can conveniently be consolidated into an effective total energy 
density, pressure, momentum density, and anisotropic stress: 


1 
Ptot = P+ JA (2 = Bnav ) + pe; (100) 
1 2 Hv H PE 
Prot = P+ Fy (2p? + 4pp + Tuun” — 4q") + = (101) 
fe} 1 V 
Qe = du + zy (Plu — Bt uvg”) + Gis (102) 
fe} 1 a 
Li = Ty + aa [-(p + 3p) Tuy — ITaly Tv + Sdent] + Ts (103) 
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These general expressions simplify in the case of a perfect fluid (or minimally coupled scalar field, 
or isotropic one-particle distribution function), i.e., for qu =0 = Tyr, to 


p pe 
asalti AE 104 
Prot pir e+e (104) 
p pe 
sapi 0o bE, 105 
Ptot ptt) + (105) 
t = qf, (106) 
tot = nE. (107) 


Note that nonlocal bulk effects can contribute to effective imperfect fluid terms even when the 
matter on the brane has perfect fluid form: There is in general an effective momentum density and 
anisotropic stress induced on the brane by massive KK modes of the 5D graviton. 

The effective total equation of state and sound speed follow from Equations and as 


_ Prot _ Wt (1+ 2w)p/2A + pe/3p 


tior 2 22 108 
tot Pot 1+ p/2à + pe/p — 
n —1 

2 _ Ptot 2, PHP Ape | | 4pe 
n E AEN T A es e E 109 
tot = ot at Nl aara me 


where w = p/p and c? = p/p. At very high energies, i.e., p > A, we can generally neglect pe (e.g., 
in an inflating cosmology), and the effective equation of state and sound speed are stiffened: 


Wtot S 2w +1, 2, eetwti. (110) 


This can have important consequences in the early universe and during gravitational collapse. For 
example, in a very high-energy radiation era, w = 1/3, the effective cosmological equation of state 
is ultra-stiff: wiot ~ 5/3. In late-stage gravitational collapse of pressureless matter, w = 0, the 
effective equation of state is stiff, Wtot + 1, and the effective pressure is nonzero and dynamically 
important. 


3.4 Conservation equations 


Conservation of Ty gives the standard general relativity energy and momentum conservation 
equations, in the general, nonlinear case: 


p+ O(p+p) + VW qu + 2A"qu +0” Tuv = 0, (111) 


. 4 = 5v v v v a 
Aiu + zO +Vupt (e+ p)Ay + V’ tu, +A” tay + Oud” — Epvaw”g* = 0. (112) 


In these equations, an overdot denotes u’”V,, © = Vu, is the volume expansion rate of the 
q , ; u a 
u” worldlines, A, = u, = Ay, is their 4-acceleration, o,, = V ,u,) is their shear rate, and 
Wu = -4 curlu, = 4 ,,) is their vorticity rate. 
On a Friedmann brane, we get 


Ay = Wp = Ow = 9, 0 = 3H, (113) 
where H = å/a is the Hubble rate. The covariant spatial curl is given by 


curl Vp = Cpap VeV, curl Wy = Eag VWP p), (114) 
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where €,ag is the projection orthogonal to u” of the 4D brane alternating tensor, and Ve is the 
projected part of the brane covariant derivative, defined by 


Vye gS PO Se es Va we: (115) 


In a local inertial frame at a point on the brane, with u” = 69, we have: 0 = Ap = wo = doy 
E0ag = curl Vo = curl Wop, and 


Vile = 6,,°5%;...6g°ViF Ft (local inertial frame), (116) 


where i, j, k = 1, 2,3. 

The absence of bulk source terms in the conservation equations is a consequence of having As 
as the only 5D source in the bulk. For example, if there is a bulk scalar field, then there is energy- 
momentum exchange between the brane and bulk (in addition to the gravitational interaction) [313] 
z3. 

Equation may be called the “nonlocal conservation equation”. Projecting along u” gives 
the nonlocal energy conservation equation, which is a propagation equation for pg. In the general, 
nonlinear case, this gives 


. 4 > v 
pe + 3Ope + V"qu +2A"q; + oN Tiu = 
1 
2A [30 at + 3(p +p)” ty, + O (2q qu + 1’ Tuv) + 6A Tuv 
-QV up + 3qhV" Tuv + 308’ V qu + 30" Tautu” — 30” qugv|. (117) 


Projecting into the comoving rest space gives the nonlocal momentum conservation equation, which 
is a propagation equation for i 


G 4 le 4 Sy HE P-E Vine va E 
GH) + 30%: + 3 Vue ot gre An +V Taw +A Taw + On Wy TEn Wla = 


1 = = Vs VS 
D [a0 + p)V pp + 6(p + DV’ Tuv +6” iiun) + 2T” Vup 


—6n°8 (Vas = Targu) = 3na V ane? + 12q’V pq — 69” vay — 244V" Qu 
—6T yao ag + 60 pot? qg + T, — bepagw mn?” gy 


+4(p + p)Oqu + 2q A” qu + 6Apg qn + 4q Tag]. (118) 


The 1 + 3-covariant decomposition shows two key features: 


e Inhomogeneous and anisotropic effects from the 4D matter-radiation distribution on the 
brane are a source for the 5D Weyl tensor, which nonlocally “backreacts” on the brane via 
its projection E,,. 


e There are evolution equations for the dark radiative (nonlocal, Weyl) energy (pe) and mo- 
mentum (Cia) densities (carrying scalar and vector modes from bulk gravitons), but there is 
no evolution equation for the dark radiative anisotropic stress (a) (carrying tensor, as well 


as scalar and vector, modes), which arises in both evolution equations. 


In particular cases, the Weyl anisotropic stress nE 


jv may drop out of the nonlocal conservation 
equations, i.e., when we can neglect ot” nE 


uvi Vrai, and AYT. This is the case when we consider 
linearized perturbations about an FRW background (which remove the first and last of these terms) 
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and further when we can neglect gradient terms on large scales (which removes the second term). 
This case is discussed in Section [6] But in general, and especially in astrophysical contexts, the 
Te terms cannot be neglected. Even when we can neglect these terms, it arises in the field 
equations on the brane. 

All of the matter source terms on the right of these two equations, except for the first term on 
the right of Equation (£18), are imperfect fluid terms, and most of these terms are quadratic in 
the imperfect quantities q, and muy. For a single perfect fluid or scalar field, only the Vip term 
on the right of Equation survives, but in realistic cosmological and astrophysical models, 
further terms will survive. For example, terms linear in 7, will carry the photon quadrupole in 
cosmology or the shear viscous stress in stellar models. If there are two fluids (even if both fluids 
are perfect), then there will be a relative velocity v, generating a momentum density qu = pup, 
which will serve to source nonlocal effects. 

In general, the 4 independent equations in Equations (117) and (118) constrain 4 of the 9 
independent components of E y on the brane. What is missing is an evolution equation for Te 
which has up to 5 independent components. These 5 degrees of freedom correspond to the 5 
polarizations of the 5D graviton. Thus in general, the projection of the 5-dimensional field equations 
onto the brane does not lead to a closed system, as expected, since there are bulk degrees of freedom 
whose impact on the brane cannot be predicted by brane observers. The KK anisotropic stress 

E 


Ty encodes the nonlocality. 


In special cases the missing equation does not matter. For example, if TE, = 0 by symmetry, as 


in the case of an FRW brane, then the evolution of €,,,, is determined by Equations and (L18). 
If the brane is stationary (with Killing vector parallel to u”), then evolution equations are not 
needed for €,,, although in general TE, will still be undetermined. However, small perturbations 
of these special cases will immediately restore the problem of missing information. 

If the matter on the brane has a perfect-fluid or scalar-field energy-momentum tensor, the local 
conservation equations and reduce to 


p+ O(p+p)=0, (119) 
Yup + (o +p)A, =0, (120) 


while the nonlocal conservation equations (117) and (118) reduce to 


4 = 
pe + 3Ope + Vig, + 2AMG, +o!’ Tiu = 0, (121) 


vi 4 wi v v va p+tp)s 
Vape + greAn +V nE, +A TE, +a, E — e wE = Ban T (122) 


1 
3 
Equation (122) shows that [391] 


e if €,, = 0 and the brane energy-momentum tensor has perfect fluid form, then the density p 


4 
sE E 
Up) + FOG, + 


must be homogeneous, Vap = 0; 


e the converse does not hold, i.e., homogeneous density does not in general imply vanishing 
eu. 


A simple example of the latter point is the FRW case: Equation (122) is trivially satisfied, 
while Equation (121) becomes 
pe +4Hpe = 0. (123) 


This equation has the dark radiation solution 


a 4 
pe = peo (2) ; (124) 
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If Euv = 0, then the field equations on the brane form a closed system. Thus for perfect fluid 
branes with homogeneous density and €,,, = 0, the brane field equations form a consistent closed 
system. However, this is unstable to perturbations, and there is also no guarantee that the resulting 
brane metric can be embedded in a regular bulk. 

It also follows as a corollary that inhomogeneous density requires nonzero Env: 


Vup4#0 > Ew#0. (125) 


For example, stellar solutions on the brane necessarily have Epy # 0 in the stellar interior if it is 
non-uniform. Perturbed FRW models on the brane also must have E v 4 0. Thus a nonzero Env, 


E, is inevitable in realistic astrophysical and cosmological models. 


and in particular a nonzero Tiy, 


3.5 Propagation and constraint equations on the brane 


The propagation equations for the local and nonlocal energy density and momentum density are 
supplemented by further 1 + 3-covariant propagation and constraint equations for the kinematic 
quantities O, A}, Wu, Cuv, and for the free gravitational field on the brane. The kinematic quantities 
govern the relative motion of neighbouring fundamental world-lines. The free gravitational field 
on the brane is given by the brane Weyl tensor C,,,93. This splits into the gravito-electric and 
gravito-magnetic fields on the brane: 


a 1 [ea 
Biv =C epg” = Bigs Hig = 5 Ena A yu" = Hwy (126) 


where E,,, is not to be confused with €,,. The Ricci identity for u” 


1 
V [uV ijua = z Rovnou’, (127) 
and the Bianchi identities 
1 
V’ Cavab = Viu (-Rae + Rava) ; (128) 


produce the fundamental evolution and constraint equations governing the above covariant quan- 
tities. The field equations are incorporated via the algebraic replacement of the Ricci tensor R,, 
by the effective total energy-momentum tensor, according to Equation (69). The brane equations 
are derived directly from the standard general relativity versions by simply replacing the energy- 
momentum tensor terms p,... by pPtot,.... For a general fluid source, the equations are given 
in [305] [229]. In the case of a single perfect fluid or minimally-coupled scalar field, the equations 
reduce to the following nonlinear equations: 


e Generalized Raychaudhuri equation (expansion propagation): 


PORR i - 2 2 
0+ 50? +0vo”” —2wpwt — VA AA+ (p +3p)—A = = (2p 3p) Í k?pe. (129) 
e Vorticity propagation: 
. 2 1 v 
Wu + 3 Own +3 curl Ay — Cuv” = 0. (130) 
e Shear propagation: 
. 2 > a k? £ 
F (uv) + zew + Buy — V Avy + Caur + Ww) = AA = > Tw: (131) 
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Gravito-electric propagation (Maxwell—Wey1 E-dot equation): 


N 


K 


Eu) tO EL yy —curl Hpv ze + DP) Ow 2A captu? — 30 0p Er)? Hw Eag En = 
Ke p 
=a (0 + P) >on 
kÊ 9 = a a 
[4peou + 390) + Om’, + 3V (ugk) + 6A Wage) + 30° Teg + 3w Eau Th | (132) 


Gravito-magnetic propagation (Maxwell-Weyl H-dot equation): 


Hy) + OA + curl yy — 30au H| + weap Hn + 2A eap En? = 


7 [curl Ti — 30 G5) + Talun Te . (133) 
Vorticity constraint: 
Vi wy, — A® wy, =0. (134) 
Shear constraint: 
We — curlw, — *¥,0 + 2Epyaw” A® = —K ge. (135) 
Gravito-magnetic constraint: 
curlopy + V (utr — Huv + 2Ay,w,) = 0. (136) 


Gravito-electric divergence (Maxwell—Wey] div-E equation): 

V” Env VN pp — Gad pH + Sg” = 
P p+ = (upe -204 — 3Ẹ°nE, + 30,748 — 9e,” wg (137) 
` uP 6 pPE du Hv H V H vla s 


Gravito-magnetic divergence (Maxwell-Weyl div-H equation): 


V” Hpv — 6 (pF wily + Eura” gE — 3B a” = 
2 
k 


5 (8pew, — 3 curl aë — Bep” *0 rEg — 37E”) . (138) 


K(p + p) Su, + 


Gauss—Codazzi equations on the brane (with w, = 0): 


Ri) + Ow) + Qouy — V wr) — Au Ar) = KTE, (139) 
9 2 
Rt + 50? — ouo” — 2n®p— 2A = ee + 2K pe, (140) 
where Riv is the Ricci tensor for 3-surfaces orthogonal to u” on the brane, and R+ = ae ore 


The standard 4D general relativity results are regained when A~! — 0 and €,,, = 0, which sets 


all right hand sides to zero in Equations (129) [30 [310 [132] (133) (134) (135) (136) [137] [138] [39 


T40). 


Together with Equations (019) [20] [21 022), these equations govern the dynamics of the 


matter and gravitational fields on the brane, incorporating both the local, high-energy (quadratic 
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Figure 4: The evolution of the dimensionless shear parameter Oshear = 07/6H? on a Bianchi I 
brane, fora V = 4m? ¢? model. The early and late-time expansion of the universe is isotropic, but 
the shear dominates during an intermediate anisotropic stage. (Figure taken from [309].) 


energy-momentum) and nonlocal, KK (projected 5D Weyl) effects from the bulk. High-energy 
terms are proportional to p/A, and are significant only when p > à. The KK terms contain pe, Lie 
and Ro with the latter two quantities introducing imperfect fluid effects, even when the matter 
has perfect fluid form. 

Bulk effects give rise to important new driving and source terms in the propagation and con- 
straint equations. The vorticity propagation and constraint, and the gravito-magnetic constraint 
have no direct bulk effects, but all other equations do. High-energy and KK energy density terms 
are driving terms in the propagation of the expansion ©. The spatial gradients of these terms 
provide sources for the gravito-electric field E „v. The KK anisotropic stress is a driving term in 
the propagation of shear o,,, and the gravito-electric/gravito-magnetic fields, E and H,,, respec- 
tively, and the KK momentum density is a source for shear and the gravito-magnetic field. The 4D 
Maxwell—Wey1 equations show in detail the contribution to the 4D gravito-electromagnetic field 
on the brane, i.e., (Euv, Huv), from the 5D Wey! field in the bulk. 

An interesting example of how high-energy effects can modify general relativistic dynamics 
arises in the analysis of isotropization of Bianchi spacetimes. For a Binachi type I brane, Equa- 
tion becomes [309 

ay (eee ee 141 
= Fe (lta) +a ee 
if we neglect the dark radiation, where a and H are the average scale factor and expansion rate, 
and X is the shear constant. In general relativity, the shear term dominates as a > 0, but in the 
brane-world, the high-energy p° term will dominate if w > 0, so that the matter-dominated early 
universe is isotropic [309] [68] [67] [416] (376) [23] [93]. This is illustrated in Figure Ø] 

Note that this conclusion is sensitive to the assumption that pe ~ 0, which by Equation 

implies the restriction 
ot ne 0. (142) 
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Relaxing this assumption can lead to non-isotropizing solutions [66]. 

The system of propagation and constraint equations, i.e., Equations (19, [20], [0210 022) 
and (129, [30] [31 [32 [733] [34 [735] (736) [37 038] [39] L40), is exact and nonlinear, applicable 
to both cosmological and astrophysical modelling, including strong-gravity effects. In general the 
system of equations is not closed: There is no evolution equation for the KK anisotropic stress 


E 
Tiv: 
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4 Gravitational Collapse and Black Holes on the Brane 


The physics of brane-world compact objects and gravitational collapse is complicated by a number 
of factors, especially the confinement of matter to the brane, while the gravitational field can access 
the extra dimension, and the nonlocal (from the brane viewpoint) gravitational interaction between 
the brane and the bulk. Extra-dimensional effects mean that the 4D matching conditions on the 
brane, i.e., continuity of the induced metric and extrinsic curvature across the 2-surface boundary, 
are much more complicated to implement [160]. High-energy corrections increase 
the effective density and pressure of stellar and collapsing matter. In particular this means that 
the effective pressure does not in general vanish at the boundary 2-surface, changing the nature 
of the 4D matching conditions on the brane. The nonlocal KK effects further complicate the 
matching problem on the brane, since they in general contribute to the effective radial pressure 
at the boundary 2-surface. Gravitational collapse inevitably produces energies high enough, i.e., 
p > À, to make these corrections significant. 

We expect that extra-dimensional effects will be negligible outside the high-energy, short- 
range regime. The corrections to the weak-field potential, Equation (41), are at the second 
post-Newtonian (2PN) level [211]. However, modifications to Hawking radiation may bring 
significant corrections even for solar-sized black holes, as discussed below. 

A vacuum on the brane, outside a star or black hole, satisfies the brane field equations 


Ruw = —Ew, Re, =0=E&*,, V" Ew = 0. (143) 


The Weyl term €,,, will carry an imprint of high-energy effects that source KK modes (as discussed 
above). This means that high-energy stars and the process of gravitational collapse will in general 
lead to deviations from the 4D general relativity problem. The weak-field limit for a static spherical 
source, Equation (41), shows that €,,, must be nonzero, since this is the term responsible for the 
corrections to the Newtonian potential. 


4.1 The black string 


The projected Weyl term vanishes in the simplest candidate for a black hole solution. This is 
obtained by assuming the exact Schwarzschild form for the induced brane metric and “stacking” 
it into the extra dimension {76}, 


Gls? = e?it dade” + dy?, (144) 
dr? 


Gav = e2ltl/2g y = —(1 — 2GM/r)dt? + ——__— 
Gu = Plg ( /r)dt" + 1—2GM/r 


+r7dn?. (145) 
(Note that Equation is in fact a solution of the 5D field equations if Juv is any 4D Einstein 
vacuum solution, i.e., if Ruy = 0, and this can be generalized to the case Ruy = —AGur [20].) 

Each {y = const.} surface is a 4D Schwarzschild spacetime, and there is a line singularity along 
r = 0 for all y. This solution is known as the Schwarzschild black string, which is clearly not 
localized on the brane y = 0. Although “)Cagcp Æ 0, the projection of the bulk Weyl tensor 
along the brane is zero, since there is no correction to the 4D gravitational potential: 


_ GM 


r 


V(r) 


The violation of the perturbative corrections to the potential signals some kind of non-AdS; pathol- 
ogy in the bulk. Indeed, the 5D curvature is unbounded at the Cauchy horizon, as y —> oo [76]: 


Env = 0. (146) 
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Furthermore, the black string is unstable to large-scale perturbations [183]. 

Thus the “obvious” approach to finding a brane black hole fails. An alternative approach is 
to seek solutions of the brane field equations with nonzero €,,, [105]. Brane solutions of static 
black hole exteriors with 5D corrections to the Schwarzschild metric have been found 
[224], but the bulk metric for these solutions has not been found. Numerical 
integration into the bulk, starting from static black hole solutions on the brane, is plagued with 
difficulties [78]. 


4.2 Taylor expansion into the bulk 


One can use a Taylor expansion, as in Equation (88), in order to probe properties of a static black 
hole on the brane [106]. (An alternative expansion scheme is discussed in [74].) For a vacuum 
brane metric, 


a ~ 2 
uv (X,Y) = Jur (x, 0) — Ep (a, 0+)y? E Zen (2 0+)|y]" 
il 32 aß m 4 
tiz Env — enw + 2RyavpE + 6E, Ear yY SP rea (148) 


y=0+ 


This shows in particular that the propagating effect of 5D gravity arises only at the fourth order 
of the expansion. For a static spherical metric on the brane, 


Judr” de” = —F (r)dt? + in + 7r7dn?, (149) 
the projected Weyl term on the brane is given by 
fm = = [r - =E], (150) 
ina p=], (151) 
Ey =-1+H+5H (+2). (152) 


These components allow one to evaluate the metric coefficients in Equation (148). For example, the 
area of the 5D horizon is determined by ggg; defining y(r) as the deviation from a Schwarzschild 
form for H, i.e., 

H(r) =1- Pa +4(r), (153) 


where m is constant, we find 
= 2 1 2 2 1 1 1 / 1 / 4 
goo(r,y) =r° — Y 1+ sly! y+ ea ae See ie) yt... (154) 


This shows how w and its r-derivatives determine the change in area of the horizon along the extra 
dimension. For the black string w = 0, and we have gge(r, y) = r°. For a large black hole, with 
horizon scale >> £, we have from Equation that 


4ml? 
3r3 ` 


yr — (155) 


This implies that ggg is decreasing as we move off the brane, consistent with a pancake-like shape of 
the horizon. However, note that the horizon shape is tubular in Gaussian normal coordinates f 
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4.3 The “tidal charge” black hole 


The equations (143) form a system of constraints on the brane in the stationary case, including 
the static spherical case, for which 


O=0=wp =T, pe =0=4i Sc |. (156) 


The nonlocal conservation equations V”€,,, = 0 reduce to 


z 4 7 
zVupe + speAy + Vi ain +A” Tiy = 0, (157) 
where, by symmetry, 
1 
iy = Ile (Fu = raro) ; (158) 


for some I¢(r), with r, being the unit radial vector. The solution of the brane field equations 
requires the input of €,, from the 5D solution. In the absence of a 5D solution, one can make an 
assumption about Epy or guv to close the 4D equations. 

If we assume a metric on the brane of Schwarzschild-like form, i.e., H = F in Equation (149), 
then the general solution of the brane field equations is [105] 


2GM  2Gé 
_2GM , 2G0Q 


F=1 An 


1 

r r ven 
2GL 

Ew =— “a [Up Ur — 2ryro + A ; (160) 


where Q is a constant. It follows that the KK energy density and anisotropic stress scalar (defined 
via Equation (158)) are given by 
LQ 1 
aero tha pile. (161) 
The solution (159) has the form of the general relativity Reissner—-Nordstrém solution, but 
there is no electric field on the brane. Instead, the nonlocal Coulomb effects imprinted by the bulk 
Weyl tensor have induced a “tidal” charge parameter Q, where Q = Q(M), since M is the source 
of the bulk Wey] field. We can think of the gravitational field of M being “reflected back” on the 
brane by the negative bulk cosmological constant [104]. If we impose the small-scale perturbative 
limit (r < £) in Equation (40), we find that 


PE 


Q=-2M. (162) 


Negative Q is in accord with the intuitive idea that the tidal charge strengthens the gravitational 
field, since it arises from the source mass M on the brane. By contrast, in the Reissner—Nordstr6m 
solution of general relativity, Q œ +q?, where q is the electric charge, and this weakens the 
gravitational field. Negative tidal charge also preserves the spacelike nature of the singularity, and 
it means that there is only one horizon on the brane, outside the Schwarzschild horizon: 


WQ 


Th = GM GM2 


1+ 4/1 =GM 


iai (163) 


The tidal-charge black hole metric does not satisfy the far-field r7’ correction to the gravi- 
tational potential, as in Equation (41), and therefore cannot describe the end-state of collapse. 
However, Equation (159) shows the correct 5D behaviour of the potential (x r~?) at short dis- 
tances, so that the tidal-charge metric could be a good approximation in the strong-field regime 
for small black holes. 


28 


4.4 Realistic black holes 


Thus a simple brane-based approach, while giving useful insights, does not lead to a realistic black 
hole solution. There is no known solution representing a realistic black hole localized on the brane, 
which is stable and without naked singularity. This remains a key open question of nonlinear 
brane-world gravity. (Note that an exact solution is known for a black hole on a 1 + 2-brane in a 
4D bulk [142], but this is a very special case.) Given the nonlocal nature of E,,,,, it is possible that 
the process of gravitational collapse itself leaves a signature in the black hole end-state, in contrast 
with general relativity and its no-hair theorems. There are contradictory indications about the 
nature of the realistic black hole solution on the brane: 


e Numerical simulations of highly relativistic static stars on the brane [430] indicate that 
general relativity remains a good approximation. 


e Exact analysis of Oppenheimer—Snyder collapse on the brane shows that the exterior is non- 
static [160], and this is extended to general collapse by arguments based on a generalized 
AdS/CFT correspondence [409] [40]. 


The first result suggests that static black holes could exist as limits of increasingly compact static 

stars, but the second result and conjecture suggest otherwise. This remains an open question. 

More recent numerical evidence is also not conclusive, and it introduces further possible subtleties 

to do with the size of the black hole 433]. 

On very small scales relative to the AdSs curvature scale, r « @, the gravitational potential 
becomes 5D, as shown in Equation (40), 

GLM G5M 

V(r) x — = —. 


r 


(164) 
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In this regime, the black hole is so small that it does not “see” the brane, so that it is approximately 
a 5D Schwarzschild (static) solution. However, this is always an approximation because of the self- 
gravity of the brane (the situation is different in ADD-type brane-worlds where there is no brane 
tension). As the black hole size increases, the approximation breaks down. Nevertheless, one might 
expect that static solutions exist on sufficiently small scales. Numerical investigations appear to 
confirm this [264] [433]: Static metrics satisfying the asymptotic AdS; boundary conditions are 
found if the horizon is small compared to £, but no numerical convergence can be achieved close 
to £. The numerical instability that sets in may mask the fact that even the very small black holes 
are not strictly static. Or it may be that there is a transition from static to non-static behaviour. 
Or it may be that static black holes do exist on all scales. 

The 4D Schwarzschild metric cannot describe the final state of collapse, since it cannot incor- 
porate the 5D behaviour of the gravitational potential in the strong-field regime (the metric is 
incompatible with massive KK modes). A non-perturbative exterior solution should have nonzero 
Euv in order to be compatible with massive KK modes in the strong-field regime. In the end-state 
of collapse, we expect an Epy which goes to zero at large distances, recovering the Schwarzschild 
weak-field limit, but which grows at short range. Furthermore, €,, may carry a Weyl “fossil 
record” of the collapse process. 


4.5 Oppenheimer—Snyder collapse gives a non-static black hole 


The simplest scenario in which to analyze gravitational collapse is the Oppenheimer—Snyder model, 
i.e., collapsing homogeneous and isotropic dust [160]. The collapse region on the brane has an FRW 
metric, while the exterior vacuum has an unknown metric. In 4D general relativity, the exterior is 
a Schwarzschild spacetime; the dynamics of collapse leaves no imprint on the exterior. 
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The collapse region has the metric 


a(r)? [dr? + r?dQ?] 


2 2 
= 1 
ds dt” + ior (165) 
where the scale factor satisfies the modified Friedmann equation (see below), 
a? 8G p pe 
ee Fall Se ea Ey 1 
z2 Sofi e+e) (166) 


The dust matter and the dark radiation evolve as 


3 4 
a a 
P = Po (2) , pe = peo (2) i (167) 
H H 


where ao is the epoch when the cloud started to collapse. The proper radius from the centre of 

the cloud is R(T) = ra(r)/(1 + 4kr?°). The collapsing boundary surface X is given in the interior 

comoving coordinates as a free-fall surface, i.e. r = ro = const., so that Rs(7) = roa(T)/(1+ krô). 
We can rewrite the modified Friedmann equation on the interior side of X as 


_2GM_ 3GM? Q 


52 
á R + A4TAR4 i R? 


FE, (168) 


where the “physical mass” M (total energy per proper star volume), the total “tidal charge” Q, 
and the “energy” per unit mass E are given by 
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M= ee, (169) 
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Now we assume that the exterior is static, and satisfies the standard 4D junction conditions. 
Then we check whether this exterior is physical by imposing the modified Einstein equations (143). 
We will find a contradiction. 

The standard 4D Darmois—Israel matching conditions, which we assume hold on the brane, 
require that the metric and the extrinsic curvature of © be continuous (there are no intrinsic 
stresses on X). The extrinsic curvature is continuous if the metric is continuous and if R is 
continuous. We therefore need to match the metrics and R across X. 

The most general static spherical metric that could match the interior metric on © is 


dR? 


2Gm(R) 
7 1—2Gm(R)/R 


R?dQ?. 172 
R + (172) 


ds? = —F(R)? i | dt? + 
We need two conditions to determine the functions F(R) and m(R). Now È is a freely falling 
surface in both metrics, and the radial geodesic equation for the exterior metric gives R? = —1 + 
2Gm(R)/R + E/F(R)?, where E is a constant and the dot denotes a proper time derivative, as 
above. Comparing this with Equation gives one condition. The second condition is easier to 
derive if we change to null coordinates. The exterior static metric, with 


dR 


dv = dt + — 
v= t Gm 
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2 
ds? = —F? (1 = =) dv? +2F dv dR + R? d?. (173) 


The interior Robertson—Walker metric takes the form [160] 


72, [L — (k + å?) R?/a?] dv? 27.,dudR 
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ds? = + R?d0?, (174) 


where 
dR 


1 
dt =T»d 1+ <—kr? | ———__. 
T= T, vt (14 gh) ae 


Comparing Equations (173) and (74) on © gives the second condition. The two conditions together 
imply that F is a constant, which we can take as F(R) = 1 without loss of generality (choosing 
E = E +1), and that 
3M? Q 
RSM tA yp, 
m(R)= M+ a CR 
In the limit \~' = 0 = Q, we recover the 4D Schwarzschild solution. In the general brane-world 
case, Equations (72) and imply that the brane Ricci scalar is 


(175) 


9GM? 
K, = —. 176 
MIRE oe) 
However, this contradicts the field equations (143), which require 
RP, = 0. (177) 


It follows that a static exterior is only possible if M/A = 0, which is the 4D general relativity limit. 
In the brane-world, collapsing homogeneous and isotropic dust leads to a non-static exterior. Note 
that this no-go result does not require any assumptions on the nature of the bulk spacetime, which 
remains to be determined. 

Although the exterior metric is not determined (see [179] for a toy model), we know that its 
non-static nature arises from 


e 5D bulk graviton stresses, which transmit effects nonlocally from the interior to the exterior, 
and 


e the non-vanishing of the effective pressure at the boundary, which means that dynamical 
information from the interior can be conveyed outside via the 4D matching conditions. 


The result suggests that gravitational collapse on the brane may leave a signature in the exterior, 
dependent upon the dynamics of collapse, so that astrophysical black holes on the brane may in 
principle have KK “hair”. It is possible that the non-static exterior will be transient, and will tend 
to a static geometry at late times, close to Schwarzschild at large distances. 


4.6 AdS/CFT and black holes on 1-brane RS-type models 


Oppenheimer-Snyder collapse is very special; in particular, it is homogeneous. One could argue 
that the non-static exterior arises because of the special nature of this model. However, the un- 
derlying reasons for non-static behaviour are not special to this model; on the contrary, the role 
of high-energy corrections and KK stresses will if anything be enhanced in a general, inhomoge- 
neous collapse. There is in fact independent heuristic support for this possibility, arising from the 
AdS/CFT correspondence. 
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The basic idea of the correspondence is that the classical dynamics of the AdSs5 gravitational 
field correspond to the quantum dynamics of a 4D conformal field theory on the brane. This 
correspondence holds at linear perturbative order [130], so that the RS 1-brane infinite AdS5 
brane-world (without matter fields on the brane) is equivalently described by 4D general relativity 
coupled to conformal fields, 

Cy = enero. (178) 


According to a conjecture [409], the correspondence holds also in the case where there is strong 
gravity on the brane, so that the classical dynamics of the bulk gravitational field of the brane 
black hole are equivalent to the dynamics of a quantum-corrected 4D black hole (in the dual 
CFT-plus-gravity description). In other words [140]: 


e Quantum backreaction due to Hawking radiation in the 4D picture is described as classical 
dynamics in the 5D picture. 


e The black hole evaporates as a classical process in the 5D picture, and there is thus no 
stationary black hole solution in RS 1-brane. 


A further remarkable consequence of this conjecture is that Hawking evaporation is dramatically 
enhanced, due to the very large number of CFT modes of order (¢/¢,)?. The energy loss rate due 
to evaporation is 


M 1 
SSN ae 1 
M (aan) , oe 
where N is the number of light degrees of freedom. Using N ~ ¢?/G, this gives an evaporation 
timescale [409] 
M \° (1mm)? 
teva, i= = e . 1 
e~ (s) Ge) ” aso 


A more detailed analysis shows that this expression should be multiplied by a factor ~ 
100. Then the existence of stellar-mass black holes on long time scales places limits on the AdS5 
curvature scale that are more stringent than the table-top limit, Equation (6). The existence of 
black hole X-ray binaries implies 

£ S 107°? mm, (181) 


already an order of magnitude improvement on the table-top limit. 

One can also relate the Oppenheimer-Snyder result to these considerations. In the AdS/CFT 
picture, the non-vanishing of the Ricci scalar, Equation (176), arises from the trace of the Hawking 
CFT energy-momentum tensor, as in Equation (178). If we evaluate the Ricci scalar at the black 
hole horizon, R ~ 2GM, using À = 6M$/M>, we find 


M; g6 


RY, aii “ME (182) 
The CFT trace on the other hand is given by T ~ N Th /MÈ , so that 
M. 12 48 
8rGTC®) ~ a (183) 


Thus the Oppenheimer—Snyder result is qualitatively consistent with the AdS/CFT picture. 

Clearly the black hole solution, and the collapse process that leads to it, have a far richer 
structure in the brane-world than in general relativity, and deserve further attention. In particular, 
two further topics are of interest: 
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e Primordial black holes in 1-brane RS-type cosmology have been investigated in 
[185] [315] QI [387]. High-energy effects in the early universe (see the next Section B) can 
significantly modify the evaporation and accretion processes, leading to a prolonged survival 
of these black holes. Such black holes evade the enhanced Hawking evaporation described 
above when they are formed, because they are much smaller than £. 


e Black holes will also be produced in particle collisions at energies = Ms, possibly well below 
the Planck scale. In ADD brane-worlds, where M4+a = O(TeV) is not ruled out by current 
observations if d > 1, this raises the exciting prospect of observing black hole production 
signatures in the next-generation colliders and cosmic ray detectors (see [75] [170] [739]). 
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5 Brane-World Cosmology: Dynamics 


A 1 + 4-dimensional spacetime with spatial 4-isotropy (4D spherical/plane/hyperbolic symmetry) 
has a natural foliation into the symmetry group orbits, which are 1 + 3-dimensional surfaces with 
3-isotropy and 3-homogeneity, i.e., FRW surfaces. In particular, the AdS; bulk of the RS brane- 
world, which admits a foliation into Minkowski surfaces, also admits an FRW foliation since it 
is 4-isotropic. Indeed this feature of 1-brane RS-type cosmological brane-worlds underlies the 
importance of the AdS/CFT correspondence in brane-world cosmology 
292 (350) [1ST]. 

The generalization of AdS; that preserves 4-isotropy and solves the vacuum 5D Einstein equa- 
tion is Schwarzschild-AdSs, and this bulk therefore admits an FRW foliation. It follows 
that an FRW brane-world, the cosmological generalization of the RS brane-world, is a part of 
Schwarzschild—AdS;, with the Z2-symmetric FRW brane at the boundary. (Note that FRW 
branes can also be embedded in non-vacuum generalizations, e.g., in Reissner-Nordstr6m—AdS; 
and Vaidya—AdS5.) 

In natural static coordinates, the bulk metric is 


dR? dr? 
Ods? = —F(R)dT? + —— + R? (| — + ran? 184 
s (R)d a0 (Set : (184) 
. PR m 
F(R) =K + F (185) 


where K = 0,+1 is the FRW curvature index, and m is the mass parameter of the black hole at 
R = 0 (recall that the 5D gravitational potential has R~? behaviour). The bulk black hole gives 
rise to dark radiation on the brane via its Coulomb effect. The FRW brane moves radially along 
the 5th dimension, with R = a(T), where a is the FRW scale factor, and the junction conditions 
determine the velocity via the Friedmann equation for a [45]. Thus one can interpret the 
expansion of the universe as motion of the brane through the static bulk. In the special case m = 0 
and da/dT = 0, the brane is fixed and has Minkowski geometry, i.e., the original RS 1-brane 
brane-world is recovered in different coordinates. 

The velocity of the brane is coordinate-dependent, and can be set to zero. We can use Gaussian 
normal coordinates, in which the brane is fixed but the bulk metric is not manifestly static [38|: 
d 2 


T Ka + r2d0?] + dy?. (186) 


lds? = —N?(t, y)dt? + A(t 
ds (tydi +A (ty) | 7R 


Here a(t) = A(t, 0) is the scale factor on the FRW brane at y = 0, and t may be chosen as proper 
time on the brane, so that N(t,0) = 1. In the case where there is no bulk black hole (m = 0), the 
metric functions are 


A 
yet (187) 
a(t) 
t 

A = a(t) | cosh (4) ed pg PD aie (EY, (188) 
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Again, the junction conditions determine the Friedmann equation. The extrinsic curvature at the 

brane is ETE. 
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Then, by Equation (68), 
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The field equations yield the first integral [38] 
At ag A 
(AA’)? a ae +m =0, (192) 


where m is constant. Evaluating this at the brane, using Equation (191), gives the modified 
Friedmann equation (194). 

The dark radiation carries the imprint on the brane of the bulk gravitational field. Thus we 
expect that E,» for the Friedmann brane contains bulk metric terms evaluated at the brane. In 
Gaussian normal coordinates (using the field equations to simplify), 


A" A ; 1 ; 
E = 3T t S Ej = (3e) ô. (193) 

Either form of the cosmological metric, Equation (184) or (86), may be used to show that 5D 
gravitational wave signals can take “short-cuts” through the bulk in travelling between points A 
and B on the brane [89] [212] [65]. The travel time for such a graviton signal is less than the time 
taken for a photon signal (which is stuck to the brane) from A to B. 

Instead of using the junction conditions, we can use the covariant 3D Gauss—Codazzi equa- 
tion (140) to find the modified Friedmann equation: 


Poh (ee a (194) 
=z 2A at 3 a?’ 
on using Equation (124), where 
k2 
m= Ë peodt. (195) 
The covariant Raychauhuri equation (129) yields 
: k2 p m K 
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which also follows from differentiating Equation (194) and using the energy conservation equation. 

When the bulk black hole mass vanishes, the bulk geometry reduces to AdSs, and pg = 0. 
In order to avoid a naked singularity, we assume that the black hole mass is non-negative, so 
that peo > 0. (By Equation (185), it is possible to avoid a naked singularity with negative m 
when K = —1, provided |m| < 4/4.) This additional effective relativistic degree of freedom is 
constrained by nucleosynthesis and CMB observations to be no more than ~ 5% of the radiation 


energy density [273] [24] [209] [47]: 
We | -< 06. (197) 


Prad [nuc 
The other modification to the Hubble rate is via the high-energy correction p/à. In order to 
recover the observational successes of general relativity, the high-energy regime where significant 
deviations occur must take place before nucleosynthesis, i.e., cosmological observations impose the 
lower limit 


à> (1 MeV) = M; > 104 GeV. (198) 
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This is much weaker than the limit imposed by table-top experiments, Equation (42). Since p?/X 
decays as a78 during the radiation era, it will rapidly become negligible after the end of the 
high-energy regime, p = X. 

If pe = 0 and K = 0 = A, then the exact solution of the Friedmann equations for w = p/p = 
const. is 


M, 
a = const. x [E(t + ta) 13D, i= 10 s, 199 
(t+) =F Sat (199) 
where w > —1. If pe #0 (but K = 0 = A), then the solution for the radiation era (w = 4) is 
V3 Mp 


a = const. x [t(t + t,)]/4, (200) 


= a Jax (1+ e/a) 


For t >> t, we recover from Equations (199) and (200) the standard behaviour, a œ #?/3(¢+1), 
whereas for t < t), we have the very different behaviour of the high-energy regime, 


p>à > ax tr, (201) 


When w = —1 we have p = po from the conservation equation. If K = 0 = A, we recover the 
de Sitter solution for pg = 0 and an asymptotically de Sitter solution for pg > 0: 


= £ A A Po = 
a = ag exp|Ho(t — to)], Hy =k 3 (1 + a), for pe = 0, (202) 
a=, p sinh[2Ho(t — to)| for pe > 0. (203) 
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A qualitative analysis of the Friedmann equations is given in [68] [67]. 


5.1 Brane-world inflation 


In 1-brane RS-type brane-worlds, where the bulk has only a vacuum energy, inflation on the brane 
must be driven by a 4D scalar field trapped on the brane. In more general brane-worlds, where the 
bulk contains a 5D scalar field, it is possible that the 5D field induces inflation on the brane via its 
effective projection [232] [196] [147] [371] [199] [198] [410] 429) [261] 277] [197] 46) [327] [222] [317] (150) [18]. 

More exotic possibilities arise from the interaction between two branes, including possible 
collision, which is mediated by a 5D scalar field and which can induce either inflation ora 
hot big-bang radiation era, as in the “ekpyrotic” or cyclic scenario [230}|216] [342] [406] [275) [319] 45], 
or in colliding bubble scenarios [40] [157] [158]. (See also [26] [98] [301] for colliding branes in an 
M theory approach.) Here we discuss the simplest case of a 4D scalar field ¢ with potential V(¢) 
(see for a review). 

High-energy brane-world modifications to the dynamics of inflation on the brane have been 
investigated [309]. Essentially, 
the high-energy corrections provide increased Hubble damping, since p > A implies that H is 
larger for a given energy than in 4D general relativity. This makes slow-roll inflation possible even 
for potentials that would be too steep in standard cosmology [206]. 

The field satisfies the Klein—Gordon equation 


6+3Hd+V"(d) =0. (204) 


In 4D general relativity, the condition for inflation, @ > 0, is $? < V(d), i.e., p < —4p, where 
p= i¢? +V and p= ig? —V. The modified Friedmann equation leads to a stronger condition 
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for inflation: Using Equation (194), with m = 0 = A = K, and Equation (204), we find that 


—— 
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a>0 => w<-3| (205) 


3 


where the square brackets enclose the brane correction to the general relativity result. As p/A > 0, 
the 4D result w < -4 is recovered, but for p > A, w must be more negative for inflation. In the 
very high-energy limit p/A — oo, we have w < -2, When the only matter in the universe is a 
self-interacting scalar field, the condition for inflation becomes 
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which reduces to ġ? < V when Po = i¢? +V <x. 
In the slow-roll approximation, we get 
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Hx &vl+ >], 207 
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The brane-world correction term V/A in Equation (207) serves to enhance the Hubble rate for 
a given potential energy, relative to general relativity. Thus there is enhanced Hubble ‘friction’ 
in Equation (208), and brane-world effects will reinforce slow-roll at the same potential energy. 
We can see this by defining slow-roll parameters that reduce to the standard parameters in the 
low-energy limit: 
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Self-consistency of the slow-roll approximation then requires €, |7| < 1. At low energies, V < A, 
the slow-roll parameters reduce to the standard form. However at high energies, V > A, the extra 
contribution to the Hubble expansion helps damp the rolling of the scalar field, and the new factors 
in square brackets become ~% A/V: 


4A 2X 
E X Egr 2] , N XS Ngr 2] ’ (211) 


where €gr, Her are the standard general relativity slow-roll parameters. In particular, this means 
that steep potentials which do not give inflation in general relativity, can inflate the brane-world 
at high energy and then naturally stop inflating when V drops below A. These models can be 
constrained because they typically end inflation in a kinetic-dominated regime and thus generate 
a blue spectrum of gravitational waves, which can disturb nucleosynthesis [288]. 
They also allow for the novel possibility that the inflaton could act as dark matter or quintessence 
at low energies [99] [314] [372] [349] [288] [8] [329] 290] [56] [383]. 


The number of e-folds during inflation, N = f Hdt, is, in the slow-roll approximation, 
Nw ey ERAP (212) 
~ Mja V' 2A 
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Brane-world effects at high energies increase the Hubble rate by a factor V/2A, yielding more 
inflation between any two values of @ for a given potential. Thus we can obtain a given number of 
e-folds for a smaller initial inflaton value ¢;. For V > A, Equation (212) becomes 


12873? V2 
~~ 3mMs [ wr (213) 


The key test of any modified gravity theory during inflation will be the spectrum of pertur- 
bations produced due to quantum fluctuations of the fields about their homogeneous background 
values. We will discuss brane-world cosmological perturbations in the next Section [6] In general, 
perturbations on the brane are coupled to bulk metric perturbations, and the problem is very 
complicated. However, on large scales on the brane, the density perturbations decouple from the 
bulk metric perturbations B73] [149]. For 1-brane RS-type models, there is no scalar 
zero-mode of the bulk graviton, and in the extreme slow-roll (de Sitter) limit, the massive scalar 
modes are heavy and stay in their vacuum state during inflation [149]. Thus it seems a reasonable 
approximation in slow-roll to neglect the KK effects carried by €,, when computing the density 
perturbations. 
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Figure 5: The relation between the inflaton mass m/M4 (Ms = Mp) and the brane tension 
(A/MŻ)" £ necessary to satisfy the COBE constraints. The straight line is the approximation 
used in Equation (220), which at high energies is in excellent agreement with the exact solution, 
evaluated numerically in slow-roll. (Figure taken from [310].) 


To quantify the amplitude of scalar (density) perturbations we evaluate the usual gauge- 
invariant quantity 


H 
C =R — —ôp, (214) 
È 
which reduces to the curvature perturbation R on uniform density hypersurfaces (69 = 0). This is 


conserved on large scales for purely adiabatic perturbations as a consequence of energy conserva- 
tion (independently of the field equations) [428]. The curvature perturbation on uniform density 
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hypersurfaces is given in terms of the scalar field fluctuations on spatially flat hypersurfaces 6¢ by 
ô 
=H ~ (215) 


If one makes the assumption that backreaction due to metric perturbations in the bulk can be 
neglected, the field fluctuations at Hubble crossing (k = aH) in the slow-roll limit are given 
by (6¢2) ~ (H/2m)’, a result for a massless field in de Sitter space that is also independent 
of the gravity theory [428]. For a single scalar field the perturbations are adiabatic and hence 
the curvature perturbation ¢ can be related to the density perturbations when modes re-enter the 
Hubble scale during the matter dominated era which is given by A? = 4(¢?)/25. Using the slow-roll 
equations and Equation (215), this gives 
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Thus the amplitude of scalar perturbations is increased relative to the standard result at a fixed 
value of @ for a given potential. 

A crucial assumption is that backreaction due to metric perturbations in the bulk can be 
neglected. In the extreme slow-roll limit this is necessarily correct because the coupling between 
inflaton fluctuations and metric perturbations vanishes; however, this is not necessarily the case 
when slow-roll corrections are included in the calculation. Previous work [252] [257] has shown 
that such bulk effects can be subtle and interesting (see also for other approaches). In 
particular, subhorizon inflaton fluctuations on a brane excite an infinite ladder of Kaluza—Klein 
modes of the bulk metric perturbations at first order in slow-roll parameters, and a naive slow-roll 
expansion breaks down in the high-energy regime once one takes into account the backreaction 
of the bulk metric perturbations, as confirmed by direct numerical simulations [20I]. However, 
an order-one correction to the behaviour of inflaton fluctuations on subhorizon scales does not 
necessarily imply that the amplitude of the inflaton perturbations receives corrections of order one 
on large scales; one must consistently quantise the coupled brane inflaton fluctuations and bulk 
metric perturbations. This requires a detailed analysis of the coupled brane-bulk system [70] [254]. 

It was shown that the coupling to bulk metric perturbations cannot be ignored in the equations 
of motion. Indeed, there are order-unity differences between the classical solutions without coupling 
and with slow-roll induced coupling. However, the change in the amplitude of quantum-generated 
perturbations is at next-to-leading order [254] because there is still no mixing at leading order 
between positive and negative frequencies when scales observable today crossed the horizon, so the 
Bogoliubov coefficients receive no corrections at leading order. The amplitude of perturbations 
generated is also subject to the usual slow-roll corrections on super-horizon scales. The next-order 
slow-roll corrections from bulk gravitational perturbations are calculated in and they are the 
same order as the usual Stewart—Lyth correction |407|. These results also show that the ratio of 
tensor-to-scalar perturbation amplitudes are not influenced by brane-bulk interactions at leading 
order in slow-roll. It is remarkable that the predictions from inflation theories should be so robust 
that this result holds in spite of the leading-order change to the solutions of the classical equations 
of motion. 

The scale-dependence of the perturbations is described by the spectral tilt 


_ din A? 
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where the slow-roll parameters are given in Equations (209) and (210). Because these slow-roll 
parameters are both suppressed by an extra factor A/V at high energies, we see that the spectral 


x —4de + 2n, (217) 
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index is driven towards the Harrison—Zel’dovich spectrum, ns + 1, as V/A — 00; however, as ex- 
plained below, this does not necessarily mean that the brane-world case is closer to scale-invariance 
than the general relativity case. 

As an example, consider the simplest chaotic inflation model V = img. Equation gives 
the integrated expansion from ¢; to f as 


Ne eE (oe =e). (218) 


The new high-energy term on the right leads to more inflation for a given initial inflaton value @;. 

The standard chaotic inflation scenario requires an inflaton mass m ~ 101? GeV to match the 
observed level of anisotropies in the cosmic microwave background (see below). This corresponds 
to an energy scale ~ 1016 GeV when the relevant scales left the Hubble scale during inflation, 
and also to an inflaton field value of order 3Mp. Chaotic inflation has been criticised for requiring 
super-Planckian field values, since these can lead to nonlinear quantum corrections in the potential. 

If the brane tension À is much below 10!° GeV, corresponding to Ms < 10! GeV, then the 
terms quadratic in the energy density dominate the modified Friedmann equation. In particular 
the condition for the end of inflation given in Equation (206) becomes ¢? < 2V. In the slow-roll 


approximation (using Equations (207) and (208)) ¢ ~ —M3/27¢, and this yields 


5 (Ms\" 
dena © DZ (=) M3. (219) 
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In order to estimate the value of when scales corresponding to large-angle anisotropies on the 
microwave background sky left the Hubble scale during inflation, we take Neobe ~ 55 in Equa- 
tion (218) and f = pena. The second term on the right of Equation (218) dominates, and we 


obtain 3 
165 / M; 
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Peobe © = (=) Ms. (220) 
Imposing the COBE normalization on the curvature perturbations given by Equation (216) requires 

8 2 445 
Asx M Pcobe x 2x107. (221) 
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Substituting in the value of cobe given by Equation (220) shows that in the limit of strong brane 
corrections, observations require 


m a 5 x 107° Ms, dcobe © 3 x 10? Ms. (222) 


Thus for Ms < 10!” GeV, chaotic inflation can occur for field values below the 4D Planck scale, 
Pcobe < Mp, although still above the 5D scale Ms. The relation determined by COBE constraints 
for arbitrary brane tension is shown in Figure] together with the high-energy approximation used 
above, which provides an excellent fit at low brane tension relative to M4. 

It must be emphasized that in comparing the high-energy brane-world case to the standard 
4D case, we implicitly require the same potential energy. However, precisely because of the high- 
energy effects, large-scale perturbations will be generated at different values of V than in the 
standard case, specifically at lower values of V, closer to the reheating minimum. Thus there are 
two competing effects, and it turns out that the shape of the potential determines which is the 
dominant effect [286]. For the quadratic potential, the lower location on V dominates, and the 
spectral tilt is slightly further from scale invariance than in the standard case. The same holds for 
the quartic potential. Data from WMAP and 2dF can be used to constrain inflationary models via 
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Figure 6: Constraints from WMAP data on inflation models with quadratic and quartic potentials, 
where R is the ratio of tensor to scalar amplitudes and n is the scalar spectral index. The high 
energy (H.E.) and low energy (L.E.) limits are shown, with intermediate energies in between, and 
the 1-o and 2-o contours are also shown. (Figure taken from [286].) 
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their deviation from scale invariance, and the high-energy brane-world versions of the quadratic and 
quartic potentials are thus under more pressure from data than their standard counterparts [286], 
as shown in Figure [6] 

Other perturbation modes have also been investigated: 


e High-energy inflation on the brane also generates a zero-mode (4D graviton mode) of tensor 
perturbations, and stretches it to super-Hubble scales, as will be discussed below. This zero- 
mode has the same qualitative features as in general relativity, remaining frozen at constant 
amplitude while beyond the Hubble horizon. Its amplitude is enhanced at high energies, 
although the enhancement is much less than for scalar perturbations [274]: 
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Equation means that brane-world effects suppress the large-scale tensor contribution 
to CMB anisotropies. The tensor spectral index at high energy has a smaller magnitude than 
in general relativity, 

ny, = —3e, (225) 


but remarkably the same consistency relation as in general relativity holds [206]: 


This consistency relation persists when Zə symmetry is dropped (and in a two-brane 
model with stabilized radion ). It holds only to lowest order in slow-roll, as in general 
relativity, but the reason for this [384| and the nature of the corrections are not settled. 


The massive KK modes of tensor perturbations remain in the vacuum state during slow-roll 
inflation [177]. The evolution of the super-Hubble zero mode is the same as in general 
relativity, so that high-energy brane-world effects in the early universe serve only to rescale 
the amplitude. However, when the zero mode re-enters the Hubble horizon, massive KK 
modes can be excited. 


e Vector perturbations in the bulk metric can support vector metric perturbations on the brane, 
even in the absence of matter perturbations (see the next Section [6). However, there is no 
normalizable zero mode, and the massive KK modes stay in the vacuum state during brane- 
world inflation [52]. Therefore, as in general relativity, we can neglect vector perturbations 
in inflationary cosmology. 


Brane-world effects on large-scale isocurvature perturbations in 2-field inflation have also been 


considered [I7]. Brane-world (p)reheating after inflation is discussed in p pal. 


5.2 Brane-world instanton 


The creation of an inflating brane-world can be modelled as a de Sitter instanton in a way that 
closely follows the 4D instanton, as shown in (155|. The instanton consists of two identical patches of 
AdS; joined together along a de Sitter brane (dS4) with compact spatial sections. The instanton 
describes the “birth” of both the inflating brane and the bulk spacetime, which are together 
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“created from nothing”, i.e., the point at the south pole of the de Sitter 4-sphere. The Euclidean 
AdSs metric is 


Olds? cig = dr? + 2 sinh? (r/£) lax? + sin? x dF) ; (227) 


where dgy is a 3-sphere, and r < ro. The Euclidean instanton interpolates between r = 0 
(“nothing”) and r = ro (the created universe), which is a spherical brane of radius 


Hy! = ésinh(ro/2). (228) 
After creation, the brane-world evolves according to the Lorentzian continuation, x > iHot + 7/2, 
(ds? = dr? + (€Ho)? sinh? (r/2) |-a0? + Hg” cosh? (Hot) dQ) (229) 


(see Figure[Z). 
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Figure 7: Brane-world instanton. (Figure taken from [155].) 


5.3 Models with non-empty bulk 


The single-brane cosmological model can be generalized to include stresses other than As in the 
bulk: 


e The simplest example arises from considering a charged bulk black hole, leading to the 
Reissner—Nordstrém AdS; bulk metric [22]. This has the form of Equation (184), with 


(230) 


where q is the “electric” charge parameter of the bulk black hole. The metric is a solution of 
the 5D Einstein-Maxwell equations, so that ©)T4, in Equation (60) is the energy-momentum 
tensor of a radial static 5D “electric” field. In order for the field lines to terminate on 
the boundary brane, the brane should carry a charge —q. Since the RNAdSs; metric is 4- 
isotropic, it is still possible to embed a FRW brane in it, which is moving in the coordinates 
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of Equation (184). The effect of the black hole charge on the brane arises via the junction 
conditions and leads to the modified Friedmann equation [22], 


2 
WH? =" (14 
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The field lines that terminate on the brane imprint on the brane an effective negative energy 
density —3q?/(«2a°), which redshifts like stiff matter (w = 1). The negativity of this term 
introduces the possibility that at high energies it can bring the expansion rate to zero and 
cause a turn-around or bounce (but see [205] for problems with such bounces). 


Apart from negativity, the key difference between this “dark stiff matter” and the dark 
radiation term m/a* is that the latter arises from the bulk Weyl curvature via the Ey 
tensor, while the former arises from non-vacuum stresses in the bulk via the F,,, tensor in 
Equation (69). The dark stiff matter does not arise from massive KK modes of the graviton. 


Another example is provided by the Vaidya-AdS; metric, which can be written after trans- 
forming to a new coordinate v = T + f dR/F in Equation (184), so that v = const. are null 
surfaces, and 


2 
(ds? = —F(R,v)dv? + 2dv dR + R? ag? , (232) 
1— Kr? 
R? 
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This model has a moving FRW brane in a 4-isotropic bulk (which is not static), with either a 
radiating bulk black hole (dm/dv < 0), or a radiating brane (dm/dv > 0) [77 [280] [279] [282]. 
The metric satisfies the 5D field equations (50) with a null-radiation energy-momentum 
tensor, 

Tap = wkaks, kak4 = 0, kau’ = 1, (234) 


where w x dm/dv. It follows that 
Fu = a Whips (235) 


In this case, the same effect, i.e., a varying mass parameter m, contributes to both €,,, and 
F, in the brane field equations. The modified Friedmann equation has the standard 1-brane 
RS-type form, but with a dark radiation term that no longer behaves strictly like radiation: 


2 
2 6 p m) 1 K 
=Ë, (1+ £) IEA- (236) 


By Equations and (234), we arrive at the matter conservation equations, 
VW’ Day = -buy (237) 


This shows how the brane loses (Y > 0) or gains (w < 0) energy in exchange with the bulk 
black hole. For an FRW brane, this equation reduces to 


p +3H(p +p) =—2v. (238) 
The evolution of m is governed by the 4D contracted Bianchi identity, using Equation (235): 


2 2 l 2> 
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For an FRW brane, this yields 


PE PE 3 2 , 
where PE = 3m(t)/(K a Ve 


A more complicated bulk metric arises when there is a self-interacting scalar field ® in the 
bulk [48]. In the simplest case, when there is no coupling between 
the bulk field and brane matter, this gives 


1 
OT ap = A4? g — gap |V(S) + ao eee (241) 


where ®(z,y) satisfies the 5D Klein—Gordon equation, 


(OQ —v'(®) =0. (242) 


The junction conditions on the field imply that 
Oy P(x, 0) = 0. (243) 


Then Equations (74) and (241) show that matter conservation continues to hold on the brane 
in this simple case: 


V’ Ty = 0. (244) 
From Equation (241) one finds that 
1 5 ag 
Fy = AK2 46 uw = Juv 3V(¢) + ag Pa, B ; (245) 
KÉ 2 
where 
$lx) = ®(2,0), (246) 


so that the modified Friedmann equation becomes 
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When there is coupling between brane matter and the bulk scalar field, then the Friedmann 


and conservation equations are more complicated [313] [21] B24 [276] 48]. 
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6 Brane-World Cosmology: Perturbations 


The background dynamics of brane-world cosmology are simple because the FRW symmetries 
simplify the bulk and rule out nonlocal effects. But perturbations on the brane immediately 
release the nonlocal KK modes. Then the 5D bulk perturbation equations must be solved in order 
to solve for perturbations on the brane. These 5D equations are partial differential equations for 
the 3D Fourier modes, with both initial and boundary conditions needed. 

The theory of gauge-invariant perturbations in brane-world cosmology has been extensively 
investigated and developed [305] [273] [24] [310] [99] [314] [372] [349] [288] [178] [274] [177] [334] [836] 
[91] [236] [267] [419] 260] [335] [420] 232) 268) [192] [124] [159] 235] [50] [128] [343] [388] B71] 2871 [85] 
and is qualitatively well understood. The key task is 
integration of the coupled brane-bulk perturbation equations. Special cases have been solved, 
where these equations effectively decouple [273] [24] [284] 283], and approximation schemes have 
been developed [401] [431] [390] [402] 403] 217 (364) (51) [202] [137] [823] [27] for the more general 
cases where the coupled system must be solved. Below we will also present the results of full 
numerical integration of the 5D perturbation equations in the RS case. 

From the brane viewpoint, the bulk effects, i.e., the high-energy corrections and the KK modes, 
act as source terms for the brane perturbation equations. At the same time, perturbations of mat- 
ter on the brane can generate KK modes (i.e., emit 5D gravitons into the bulk) which propagate 
in the bulk and can subsequently interact with the brane. This nonlocal interaction amongst the 
perturbations is at the core of the complexity of the problem. It can be elegantly expressed via 
integro-differential equations [336], which take the form (assuming no incoming 5D gravita- 
tional waves) 


Ax (t) = fu G(t, t) Be (t’), (248) 


where G is the bulk retarded Green’s function evaluated on the brane, and Ak, Bk are made 
up of brane metric and matter perturbations and their (brane) derivatives, and include high- 
energy corrections to the background dynamics. Solving for the bulk Green’s function, which then 
determines G, is the core of the 5D problem. 

We can isolate the KK anisotropic stress m}, as the term that must be determined from 5D 
equations. Once TE, is determined in this way, the perturbation equations on the brane form a 
closed system. The solution will be of the form (expressed in Fourier modes): 


E 


mé (t) x fæ G(t, t) Fi,(t’), (249) 


where the functional Fẹ will be determined by the covariant brane perturbation quantities and their 
derivatives. It is known in the case of a Minkowski background [377], but not in the cosmological 
case. 

The KK terms act as source terms modifying the standard general relativity perturbation 
equations, together with the high-energy corrections. For example, the linearization of the shear 
propagation equation (131) yields 

-2 7 Kp a 
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(1+ 3w)\ Tuv: (250) 
In 4D general relativity, the right hand side is zero. In the brane-world, the first source term on the 
right is the KK term, and the second term is the high-energy modification. The other modification 
is a straightforward high-energy correction of the background quantities H and p via the modified 
Friedmann equations. 

As in 4D general relativity, there are various different, but essentially equivalent, ways to 
formulate linear cosmological perturbation theory. First we describe the covariant brane-based 
approach. 
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6.1 1+ 3-covariant perturbation equations on the brane 


In the 1 + 3-covariant approach [305] [284] [307], perturbative quantities are projected vectors, 

Va = Viz), and projected symmetric tracefree tensors, Wy, = Wuv), which are gauge-invariant 

since they vanish in the background. These are decomposed into (3D) scalar, vector, and tensor 
modes as 

Va = Vu + Vu, (251) 

Wuv = ViuViyW + Su Way + Wuv, (252) 


where Wy, = Wuv) and an overbar denotes a (3D) transverse quantity, 
VeV, =0=V" Wy. (253) 


In a local inertial frame comoving with u”, i.e., u” = (1, 0), all time components may be set to 
zero: V, = (0,Vi), Won = 0, Vu = (0, Vi). 
Purely scalar perturbations are characterized by the fact that vectors and tensors are derived 
from scalar potentials, i.e., 
Va =W, =W =0. (254) 


Scalar perturbative quantities are formed from the potentials via the (3D) Laplacian, e.g., V = 
VNV aV = WV. Purely vector perturbations are characterized by 


=v. Ww = Fiu Woy, cul Vf = —2fwp, (255) 
where w,, is the vorticity, and purely tensor by 
Val =0= Vis Wuv = Wv. (256) 


The KK energy density produces a scalar mode Vape (which is present even if pe = 0 in the 
background). The KK momentum density carries scalar and vector modes, and the KK anisotropic 
stress carries scalar, vector, and tensor modes: 


E =Vid +E, (257) 
Thy = NaNe EN + Thy: (258) 


Linearizing the conservation equations for a single adiabatic fluid, and the nonlocal conservation 
equations, we obtain 


p+ O(p +p) =, (259) 
2%, + (9+ PAu =0, (260) 
. 4 = 
pe + Ope + Vig. = (261) 
. ls 4 $y +P)s 
Gi +4Hqa + 3V upe + ZpeAu +V”Ti = LEDS, 5 (262) 
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Linearizing the remaining propagation and constraint equations leads to 


~ 1 > 1 2 
Ò + 50? — VHA, + 3? + 3p) — A= -5 (2p + 3p)È — pe, (263) 
1 
Wy, + 2Hwy, + 5 curl A,, = 0, (264) 
> K2 
pv + 2H + Ep — Viu Av) = ae (265) 
; k2 k2 p 
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-5 [4peow +34E, +3HnE, +3Vi.05)|, (266) 
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Hpo +3HH,, + curl Epy = > curl ré,,, (267) 
V4w, = 0, (268) 

2, Do 
V" ou — curlw, — 3V9 = —, (269) 
curlopy + Ver) — Hpv = 0, (270) 

Z Ke Kp = Pa 7 
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Equations (259), (261), and do not provide gauge-invariant equations for perturbed quanti- 
ties, but their spatial gradients do. 

These equations are the basis for a 1 + 3-covariant analysis of cosmological perturbations from 
the brane observer’s viewpoint, following the approach developed in 4D general relativity (for a 
review, see [138]). The equations contain scalar, vector, and tensor modes, which can be separated 
out if desired. They are not a closed system of equations until Lae is determined by a 5D analysis 
of the bulk perturbations. An extension of the 1 + 3-covariant perturbation formalism to 1 + 4 
dimensions would require a decomposition of the 5D geometric quantities along a timelike extension 
u^ into the bulk of the brane 4-velocity field u”, and this remains to be done. The 1 + 3-covariant 
perturbation formalism is incomplete until such a 5D extension is performed. The metric-based 
approach does not have this drawback. 


6.2 Metric-based perturbations 


An alternative approach to brane-world cosmological perturbations is an extension of the 4D 
metric-based gauge-invariant theory 331). A review of this approach is given in 365|. In 
an arbitrary gauge, and for a flat FRW background, the perturbed metric has the form 

—2N2y A?(0;B — Si) Na 
gap = | A2(0;B—S;) A? {2Rbiz + 20;0j;C + 20GF) + fiz} | AB- x) |, (273) 


| Na A?(O;8 — xj) 2v 


where the background metric functions A, N are given by Equations ([87] 188). The scalars 
Y, R,C, a, 8,v represent scalar perturbations. The vectors S;, F;, and y; are transverse, so that 
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they represent 3D vector perturbations, and the tensor fij is transverse traceless, representing 3D 
tensor perturbations. 
In the Gaussian normal gauge, the brane coordinate-position remains fixed under perturbation, 


Ods? = [a (x,y) + guv (z, y)| dada” + dy”, (274) 


where go) is the background metric, Equation (186). In this gauge, we have 


a=B=v=\y;,=0. (275) 
In the 5D longitudinal gauge, one gets 
—-B4+C=0=-64C'. (276) 


In this gauge, and for an AdS; background, the metric perturbation quantities can all be expressed 
in terms of a “master variable” Q which obeys a wave equation 336|. In the case of scalar 
perturbations, we have for example 


1 tix A 
T (a"- 38-0), (277) 


with similar expressions for the other quantities. All of the metric perturbation quantities are 
determined once a solution is found for the wave equation 


fay As KNN N NV 
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The junction conditions (68) relate the off-brane derivatives of metric perturbations to the 
matter perturbations: 


R= 


1 1 
Əy Öguv = —K2 före +5 (à = 7) Sguv — Zaar] , (279) 
where 
6T° = —ôp, (280) 
ôT’; =a"q; (281) 
ôT’; = bpd"; + ôT*;. (282) 


For scalar perturbations in the Gaussian normal gauge, this gives 


a2 
Ayib(w,0) = 2 (25p + 36p), (283) 
ð,B(z,0) = n26p, (284) 
ae 
AyC(x,0) = on, (285) 
2 
dyR(x,0) = -2p — ð; C(x, 0), (286) 


where dz is the scalar potential for the matter anisotropic stress, 


Oni; = O,0;07 = oi 3o Orr. (287) 
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The perturbed KK energy-momentum tensor on the brane is given by 


66° = kK Ope, (288) 

6£°,; = —K*a" dé, (289) 
i K? i Ei 

ôE j= -5 SPE o] jT om j: (290) 


The evolution of the bulk metric perturbations is determined by the perturbed 5D field equa- 
tions in the vacuum bulk, 
6 G4 p =0. (291) 


Then the matter perturbations on the brane enter via the perturbed junction conditions (279). 
For example, for scalar perturbations in Gaussian normal gauge, we have 
| f (292) 


1 , 2 | ‘ 
sat, =a.) v4 (4-5) va Aa as (54-2) 0 


2N? 
For tensor perturbations (in any gauge), the only nonzero components of the perturbed Einstein 
tensor are 


; 1f 1 x y Fee DN wee N! AN on 
ô OG; = i oe ihitge(F-94) fie (Teed) ul (293) 


In the following, we will discuss various perturbation problems, using either a 1 + 3-covariant 
or a metric-based approach. 


6.3 Density perturbations on large scales 


In the covariant approach, we define matter density and expansion (velocity) perturbation scalars, 
as in 4D general relativity, 


2 
A=, Z=aV"O. (294) 
p 
Then we can define dimensionless KK perturbation scalars |805), 
a? s9 aso E Logg 
U = —V* pe, Q=-V¢, Il = —-V*r", (295) 
p p p 


where the scalar potentials qf and 7f are defined by Equations 258). The KK energy density 
(dark radiation) produces a scalar fluctuation U which is present even if pg = 0 in the background, 
and which leads to a non-adiabatic (or isocurvature) mode, even when the matter perturbations 
are assumed adiabatic |178]. We define the total effective dimensionless entropy Stot via 


Ptot Stot = a? V? Pitot = Gar" pees (296) 
where c2 = Ptot/Ptot is given in Equation (109). Then 


2_1 2 2 
9 [6 - 3 + (5 +w + ci) p/À] SPER iat a (297) 


“or = BTEC) + Ape] pllBw + 3U + 2w) F pelo) [Bp 


If pe = 0 in the background, then U is an isocurvature mode: Stot x (1+w)U. This isocurvature 
mode is suppressed during slow-roll inflation, when 1+ w = 0. 
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If pe Æ 0 in the background, then the weighted difference between U and A determines the 
isocurvature mode: Stot X (4p¢/3p)A — (1+ w)U. At very high energies, p >> A, the entropy is 
suppressed by the factor A/p. 

The density perturbation equations on the brane are derived by taking the spatial gradients of 
Equations (259), (261), and (263), and using Equations and (262). This leads to [178] 


A =3wHA — (14 w)Z, (298) 
Ż = —2HZ = (=) V2A — 2 pU — sep h + (4+ 3w)£ -— (=) — A, (299) 
Ù = (3w —1)HU + (=) (=) HA- (2) Z—aV°Q, (300) 
Ò = (3w—1)HQ =U 2g $ i G) T ee we A. (301) 


The KK anisotropic stress term II occurs only via its Laplacian, V°IL. If we can neglect this 
term on large scales, then the system of density perturbation equations closes on super-Hubble 
scales [305]. An equivalent statement applies to the large-scale curvature perturbations [273]. KK 
effects then introduce two new isocurvature modes on large scales (associated with U and Q), and 
they modify the evolution of the adiabatic modes as well [178] [284]. 

Thus on large scales the system of brane equations is closed, and we can determine the density 
perturbations without solving for the bulk metric perturbations. 

We can simplify the system as follows. The 3-Ricci tensor defined in Equation (140) leads to a 
scalar covariant curvature perturbation variable, 


C =a'V?Rt = —407HZ + 2670p (1 + £) A+ 2x70" pU. (302) 
It follows that C is locally conserved (along u” flow lines): 
C=O, M=0. (303) 


We can further simplify the system of equations via the variable 
® = K2a7pA. (304) 


This should not be confused with the Bardeen metric perturbation variable ® p, although it is the 
covariant analogue of ®y in the general relativity limit. In the brane-world, high-energy and KK 
effects mean that ®y is a complicated generalization of this expression [284] involving II, but the 
simple ® above is still useful to simplify the system of equations. Using these new variables, we 
find the closed system for large-scale perturbations: 


. Kp p a? rtp? k?p 
®=-H|1+(1 —, (14+-~)/®-]j(1 1 — 
| +( + wv) ( +4)| K +w) In Jos ( + wT l Co (305) 
2K? pe 2pe p 6c? H? pe 
U =-H]|1-3 U- 1+5- -l ®+ |] ©. 306 
| i 3H? | 3a? Hp X à (+w)k?p i 3a2Hp| ° on) 
If there is no dark radiation in the background, pg = 0, then 
=the (- J a- swan ) , (307) 
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Figure 8: The evolution of the covariant variable ®, defined in Equation (and not to be 
confused with the Bardeen potential), along a fundamental world-line. This is a mode that is well 
beyond the Hubble horizon at N = 0, about 50 e-folds before inflation ends, and remains super- 
Hubble through the radiation era. A smooth transition from inflation to radiation is modelled by 
w = $[(2— 3) tanh(N — 50) — (1 — e)], where € is a small positive parameter (chosen as € = 0.1 in 
the plot). Labels on the curves indicate the value of pọ/A, so that the general relativistic solution 


is the dashed curve (p9/A = 0). (Figure taken from [I78].) 
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and the above system reduces to a single equation for 6. At low energies, and for constant w, the 
non-decaying attractor is the general relativity solution, 


3(1 +w) 


Wea guy or (308) 


~ 
low ~ 


At very high energies, for w > —t, we get 


3 Co 2Uo 
Phish > —-—(1 — ; 309 
ne Tae Bon en) 
where Up = k2a2p9Uo, so that the isocurvature mode has an influence on ©. Initially, ® is 


suppressed by the factor A/po, but then it grows, eventually reaching the attractor value in Equa- 
tion (808). For slow-roll inflation, when 1 + w ~ e, with 0 < «<1 and H™Jé| = |e’| < 1, we get 


nigh ~ 3 A Cede, (310) 
2 po 

where N = In(a/ao), so that ® has a growing-mode in the early universe. This is different from 
general relativity, where ® is constant during slow-roll inflation. Thus more amplification of ® can 
be achieved than in general relativity, as discussed above. This is illustrated for a toy model of 
inflation-to-radiation in Figure[8] The early (growing) and late time (constant) attractor solutions 
are seen explicitly in the plots. 

The presence of dark radiation in the background introduces new features. In the radiation era 
(w = $), the non-decaying low-energy attractor becomes [188] 


Diow Z Sa —a), (311) 
a= < 0.05. (312) 
p 


The dark radiation serves to reduce the final value of ®, leaving an imprint on ®, unlike the pe = 0 
case, Equation (808). In the very high energy limit, 


r [2 As A\? | Gy. 40 
Phigh > A | Co — 0 + 16a (=) 2 0 
0 


: 1 
273 539 a) 


9 7 Po 


Thus ® is initially suppressed, then begins to grow, as in the no-dark-radiation case, eventually 
reaching an attractor which is less than the no-dark-radiation attractor. This is confirmed by the 
numerical integration shown in Figure [9] 


6.4 Curvature perturbations and the Sachs—Wolfe effect 


The curvature perturbation R on uniform density surfaces is defined in Equation (273). The 
associated gauge-invariant quantity 
dp 
¢=R+—->— (314) 
3(p +p) 
may be defined for matter on the brane. Similarly, for the Weyl “fluid” if pe 4 0 in the background, 
the curvature perturbation on hypersurfaces of uniform dark energy density is 


Ce = R+ —. (315) 
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Figure 9: The evolution of ® in the radiation era, with dark radiation present in the background. 
(Figure taken from [188].) 


On large scales, the perturbed dark energy conservation equation is 
(Spe) +4Hdpe + 4peR = 0, (316) 


which leads to 


és =0. (317) 


For adiabatic matter perturbations, by the perturbed matter energy conservation equation, 


(dp) +3H(6p+ 6p) + 3(p+p)R = 0, (318) 


we find 


¢=0. (319) 
This is independent of brane-world modifications to the field equations, since it depends on en- 
ergy conservation only. For the total, effective fluid, the curvature perturbation is defined as 
follows [273]: If pe # 0 in the background, we have 


4pe 
T |e he 320 
Soe = |B py + ol) + Ape | SE S) = 
and if pg = 0 in the background, we get 
Ope 
en a) o 321 
GAMBEI Fpl) ae 
ôC 
Ope = E (322) 


where 6C¢ is a constant. It follows that the curvature perturbations on large scales, like the density 
perturbations, can be found on the brane without solving for the bulk metric perturbations. 

Note that ¢ tot Æ 0 even for adiabatic matter perturbations; for example, if pe = 0 in the 
background, then 


1 


: o] 
Ctot =H a m 3 FE 


) (e+p) ++p/A) 


The KK effects on the brane contribute a non-adiabatic mode, although Ĉtot — 0 at low energies. 

Although the density and curvature perturbations can be found on super-Hubble scales, the 
Sachs—Wolfe effect requires ti in order to translate from density/curvature to metric perturba- 
tions. In the 4D longitudinal gauge of the metric perturbation formalism, the gauge-invariant 
curvature and metric perturbations on large scales are related by 


(323) 


Cot = R — — (3 — v) ; (324) 
R+w = —rk?a sre, (325) 


where the radiation anisotropic stress on large scales is neglected, as in general relativity, and ông is 
the scalar potential for Tees equivalent to the covariant quantity II defined in Equation (295). In 4D 
general relativity, the right hand side of Equation (325) is zero. The (non-integrated) Sachs—Wolfe 


formula has the same form as in general relativity: 


ôT 


T an (Gad TP wp = R) dec: (326) 
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The brane-world corrections to the general relativistic Sachs—Wolfe effect are then given by 


ôT (ôT 
P p T r 3 Pcdm a®/2 


rai 2 2 
) ae ( om ) Se — k’°a sne + | daa’? Ome, (327) 
g 


where Se is the KK entropy perturbation (determined by dp¢). The KK term rg cannot be 
determined by the 4D brane equations, so that 6T/T cannot be evaluated on large scales without 
solving the 5D equations. (Equation has been generalized to a 2-brane model, in which the 
radion makes a contribution to the Sachs—Wolfe effect [246].) 

The presence of the KK (Weyl, dark) component has essentially two possible effects: 


e A contribution from the KK entropy perturbation S¢ that is similar to an extra isocurvature 
contribution. 


e The KK anisotropic stress ông also contributes to the CMB anisotropies. In the absence 
of anisotropic stresses, the curvature perturbation Ġtot would be sufficient to determine the 
metric perturbation R and hence the large-angle CMB anisotropies via Equations (324] [325] 
326). However, bulk gravitons generate anisotropic stresses which, although they do not 
affect the large-scale curvature perturbation tot, can affect the relation between ¢, 4, R, and 
w, and hence can affect the CMB anisotropies at large angles. 


A simple phenomenological approximation to ông on large scales is discussed in [24], and the 
Sachs—Wolfe effect is estimated as 


ôT (=) ( A i oon (328) 
T P Jin \tdec In(teq/ta) J’ 
where t4 is the 4D Planck time, and tin is the time when the KK anisotropic stress is induced on 
the brane, which is expected to be of the order of the 5D Planck time. 
A self-consistent approximation is developed in [247], using the low-energy 2-brane approxi- 
mation [401] [431] [390] [402] [403] to find an effective 4D form for €,,,, and hence for dg. This is 


discussed below. In a single brane model in the AdS bulk, full numerical simulations were done to 
find the behaviour of d7¢ [69], as will be discussed in the next subsection. 


6.5 Full numerical solutions 


In order to study scalar perturbations fully, we need to numerically solve the coupled bulk and 
brane equations for the master variable Q. A thorough analysis was done in [69]. For this purpose, 
it is convenient to use the static coordinate where the bulk equation is simple and consider a 


moving FRW brane; 
2 


(5) qs? = E (undo! dx” + dz”). (329) 
Z 


The bulk master variable satisfies the following wave equation (see Equation (278) ) 


Pl PQl 300 1 
0 = -— += +- — =k?) Q. 330 
ðr? ðz z Oz (= ) ee) 
From the junction condition, Q satisfies a boundary condition on the brane 
1 p 6pa 
nQ2+—(14+—)Q+— A} =0. 1 
la £ ( 7 Ak? I, oon) 
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where ôn is the derivative orthogonal to the brane 
1 o o 
On = — | -Hl—+V14+ HV — }, (332) 
a OT Oz 


and A is the density perturbation in the comoving gauge. The subscript b implies that the 
quantities are evaluated on the brane. On a brane, A satisfies a wave equation 


aa dA 3pa 30a? ki(1 + w)R» 

tl 2 _ 6w)Ha— 2p? A Bia 333 

di? + (14+ 3c; — 6w)Ha i + [e + 7 + ee Tr (333a) 
A = 6c? —1— 8w + 3w’, B = 3 — 9w — 4, (333b) 


where we consider a perfect fluid on a brane with an equation state w and cs is a sound speed 
for perturbations. The above ordinal differential equation, the bulk wave equation and the 
boundary condition (831) comprise a closed set of equations for A and Qe. 

On a brane, we take the longitudinal gauge 


ds? = —(1 + 2Y)dt? + (1 + 2R)ôijdx dat. (334) 


Using the expressions for metric perturbations in terms of the master variable Q (see Eq. (277)) 
and the junction condition Eq. (831), y and R are written in terms of A and Q as 


3a°plp +A) 3H?a? + k? H d% 

k=- per O Gla? b — a? dn (Bapa) 

= 3pa?(3wp + 4p + dA (3w + 4)p? a (5 + 3w)p n k? 4 3H dù% 1 dP% 
E REX WaN? aA 3l | P T Ua? dy Ua? dn ` 
(335b) 

Other quantities of interest are the curvature perturbation on uniform density slices, 
HaV A 1 3pa?(wA-A-p)] A Ha dA k? 

= R- — + —— = |= = | 4 hh 336 

sake a Ge, É PEX ma riea ae 


where the velocity perturbation V is also written by A and Q». There are two independent 
numerical codes that can be used to solve for A and Qp. The first is the pseudo-spectral (PS) 
method used in [201] and the second is the characteristic integration (CI) algorithm developed in 


Figure [10] shows the output of the PS and CI codes for a typical simulation of a mode with 
p/X = 50 at the horizon re-enter. As expected we have excellent agreement between the two codes, 
despite the fact that they use different initial conditions. Note that for all simulations, we recover 
that A and ¢ are phase-locked plane waves, 


kn a 

A(n) x cos A(n) ~ 4¢(n); (337) 
at sufficiently late times kn >> 1, which is actually the same behaviour as seen in GR. Figure [I] 
illustrates how the ordinary superhorizon behaviour of perturbations in GR is recovered for modes 
entering the Hubble horizon in the low energy era. We see how A, y and R smoothly interpolate 
between the non-standard high-energy behaviour to the usual expectations in GR. Also shown 
in this plot is the behaviour of the KK anisotropic stress, which steadily decays throughout the 

simulation. These results confirm that at low energies, we recover GR solutions smoothly. 
At high energies p > A, there are two separate effects to consider: First, there is the modifica- 
tion of the universe’s expansion at high energies and the O(p/A) corrections to the perturbative 
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Figure 10: Comparison between typical results of the PS and CI codes for various brane quantities 
(left); and the typical behaviour of the bulk master variable (right) as calculated by the CI method. 
Very good agreement between the two different numerical schemes is seen in the left panel, despite 
the fact that they use different initial conditions. Also note that on subhorizon scales, A and 
¢ undergo simple harmonic oscillations, which is consistent with the behaviour in GR. The bulk 
profile demonstrates our choice of initial conditions: We see that the bulk master variable Q is 
essentially zero during the early stages of the simulation, and only becomes “large” when the mode 
crosses the horizon. Figure taken from [69]. 


equations of motion. Second, there is the effect of the bulk degrees of freedom encapsulated by 
the bulk master variable Q (or, equivalently, the KK fluid €,,). To separate out the two effects, 
it is useful to introduce the 4-dimensional effective theory where all O(p/X) corrections to GR 
are retained, but the bulk effects are removed by artificially setting Q = 0. In the case of radia- 
tion domination, we obtain equations for the effective theory density contrast A, and curvature 
perturbation ¢,.-¢ from Equations and with Qp = 0: 


PA k? 4pa? 18p%a? 
0 = —— 4 (= pA — Ayers (338a) 


1 3pa? 9p? a? 3HadA,.; 
= fo 4— 4 Ai : 
Coe G + yee + ET dn 


(338b) 


These give a closed set of ODEs on the brane that describe all of the O(p/A) corrections to GR. 
Since in any given model we expect the primordial value of the curvature perturbation to be 
fixed by inflation, it makes physical sense to normalize the waveforms from each theory such that 
Cso © Ger © Con © 1 for a << a. We can define a set of “enhancement factors”, which are functions 
of k that describe the relative amplitudes of A after horizon crossing in the various theories. 
Let the final amplitudes of the density perturbation with wavenumber k be Cs»(k), Cre(k) and 
Cer(k) for the 5-dimensional, effective and GR theories, respectively, given that the normalization 
Csp © Ger © Can © 1 holds. Then, we define enhancement factors as 
Cyer(k) Cso(k) 


i i Calk) 
Querlk) = Con(k)’ Qe(k) - Calk)” 7 Car(k) 


It follows that Q,..(&) represents the O(p/A) enhancement to the density perturbation, Qe (k) gives 


Qsn(k) 


(339) 
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Figure 11: The simulated behaviour of a mode on superhorizon scales. 
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On the left we show 


how the A gauge invariant switches from the high-energy behaviour predicted to the familiar GR 
result as the universe expands through the critical epoch. We also show how the KK anisotropic 
stress K267¢ steadily decays throughout the simulation, which is typical of all the cases we have 
investigated. On the right, we show the metric perturbations ~(= Y) and R(= ®) as well as the 
curvature perturbation ¢. Again, note how the GR result ® ~ —W ~ —2¢/3 is recovered at low 


energy. Figure taken from |69)}. 
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the magnification due to KK modes, while Q;5(k) gives the total 5-dimensional amplification over 
the GR case. They all increase as the scale is decreased, and that they all approach unity for k > 0. 
Since Q = 1 implies no enhancement of the density perturbations over the standard result, this 
means we recover general relativity on large scales. For all wavenumbers we see Q.e > Qe > 1, 
which implies that the amplitude magnification due to the O(p/A) corrections is always larger 
than that due to the KK modes. Interestingly, the Q-factors appear to approach asymptotically 
constant values for large k: 


OQue(k) 3.0, Qe(k)¥ 2.4, Qsv(k) #71, k®> ke, (340) 


where ke is the comoving wavenumber of the mode that enters the horizon when H = 4t. 

In cosmological perturbation theory, transfer functions are very important quantities. They 
allow one to transform the primordial spectrum of some quantity set during inflation into the 
spectrum of another quantity at a later time. In this sense, they are essentially the Fourier 
transform of the retarded Green’s function for cosmological perturbations. There are many different 
transfer functions one can define, but for our case it is useful to consider a function T(k) that will 
tell us how the initial spectrum of curvature perturbations Pe maps onto the spectrum of density 
perturbations Pa at some low energy epoch within the radiation era. It is customary to normalize 
transfer functions such that T (k; n) > 1(k > 0), which leads us to the following definition 


9 k iF A; (n) 


rb) =F laa cr ey 


Here, Ç} is the primordial value of the curvature perturbation and A(n) is the maximum ampli- 
tude of the density perturbation in the epoch of interest. As demonstrated in Figure [I] we know 
that we recover the GR result in the extreme small scale limit (k — 0), which gives the transfer 
function the correct normalization. In the righthand panel of Figure [12| we show the transfer 
functions derived from GR, the effective theory and the 5-dimensional simulations. As expected, 
the T(k;7) for each formulation match one another on subcritical scales k < ke. However, on 
supercritical scales we have Tsp > Tice > Tor- 

Note that if we are interested in the transfer function at some arbitrary epoch in the low-energy 
radiation regime Ha > ke, it is approximately given in terms of the enhancement factor as follows: 


1, k < 3Ha, 
Taolkin) = 7 ‘ (342) 
(3Ha/k)°Q;p(k), k > 3Ha, 
Now, the spectrum of density fluctuations at any point in the radiation era is given by 
Palkin = ÉT (A) Pa) (343) 
A\R;1)) = 81 31) Ha a 5 


Using Equation (842), we see that the RS matter power spectrum (evaluated in the low-energy 
regime) is ~ 50 times bigger than the GR prediction on scales given by k ~ 10%kc. 

The amplitude enhancement of perturbations is important on comoving scales < 10 AU, which 
are far too small to be relevant to present-day/cosmic microwave background measurements of 
the matter power spectrum. However, it may have an important bearing on the formation of 
compact objects such as primordial black holes and boson stars at very high energies, i.e. the 
greater gravitational force of attraction in the early universe will create more of these objects 
than in GR (different aspects of primordial black holes in RS cosmology in the context of various 


effective theories have been considered in [186] [185] 91] [3877] [386] [385]). 


60 


enhancement factors 


log(transfer function) 


© 5D simulation: log T;»(k; n) 
—$— effective theory: log Telk; n) 
general relativity: log Tor (k; n) 


-8- 


3° 2 1 0 1 2 3 


relative wavenumber log(k/ke) relative wavenumber log(k/ke) 


Figure 12: Density perturbation enhancement factors (left) and transfer functions (right) from 
simulations, effective theory, and general relativity. All of the Q factors monotonically increase 
with k/kc, and we see that the A amplitude enhancement due to O(p/A) effects Q.e is generally 
larger than the enhancement due to KK effects Qe. For asymptotically small scales k > ke, 
the enhancement seems to level off. The transfer functions in the right panel are evaluated at a 
given subcritical epoch in the radiation dominated era. The T functions show how, for a fixed 
primordial spectrum of curvature perturbations Pf, the effective theory predicts excess power 
in the A spectrum Pa œx Pps on supercritical/subhorizon scales compared to the GR result. 
The excess small-scale power is even greater when KK modes are taken into account, as shown by 
Tsp(k;7). Figure taken from [69]. 
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6.6 Vector perturbations 


The vorticity propagation equation on the brane is the same as in general relativity, 
. 1 
Wy + 2Hwy, = E curl A,,. (344) 


Taking the curl of the conservation equation (112) (for the case of a perfect fluid, qu = 0 = Tuv), 
and using the identity in Equation (255), one obtains 


curl A, = -6H fw, (345) 
as in general relativity, so that Equation (344) becomes 
üp + (2— 3c?) Hw, = 0, (346) 


which expresses the conservation of angular momentum. In general relativity, vector perturbations 
vanish when the vorticity is zero. By contrast, in brane-world cosmology, bulk KK effects can source 
vector perturbations even in the absence of vorticity [306]. This can be seen via the divergence 
equation for the magnetic part H,,, of the 4D Weyl tensor on the brane, 


oe p 4 1 = 
ZA, = 2k?(p +p) [1 + 2] wu + EK? pew, — 5K? curl gE, (347) 
where Hy, = Viu Hn. Even when w, = 0, there is a source for gravimagnetic terms on the brane 
from the KK quantity curl që. 

We define covariant dimensionless vector perturbation quantities for the vorticity and the KK 
gravi-vector term: 


āu = awp Bp =—curld. (348) 
p 
On large scales, we can find a closed system for these vector perturbations on the brane [306]: 
Gy, + (1 — 3c3) Ha, = 0, (349) 


age = 2 
By + (1 3w)HB, = SH |4 (3c — 1) a -9(1+ wÈ Gy. (350) 


Thus we can solve for @,, and Bu on super-Hubble scales, as for density perturbations. Vorticity in 
the brane matter is a source for the KK vector perturbation 8, on large scales. Vorticity decays 
unless the matter is ultra-relativistic or stiffer (w > 3), and this source term typically provides a 


decaying mode. There is another pure KK mode, independent of vorticity, but this mode decays 
like vorticity. For w = p/p = const., the solutions are 


a \ 387! 
Qu = On (=). (351) 
E a \ 8Y- TE 2(3w—1) oTa —4 

Bu = C (=) + by, (=) +2(1 + ew (=) | ; (352) 


3p0 \ ao 
where b, = 0 = é. 
Inflation will redshift away the vorticity and the KK mode. Indeed, the massive KK vector 
modes are not excited during slow-roll inflation [53] [365]. 
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6.7 Tensor perturbations 


The covariant description of tensor modes on the brane is via the shear, which satisfies the wave 
equation [306 


= 


» . 1 ; 
Vow — Cw — 5H pv [2a 5 (p 3p (0+ 3p)£) Ou — K? (hE, +2H7E,). (353) 


Unlike the density and vector perturbations, there is no closed system on the brane for large scales. 
The KK anisotropic stress Tiv is an unavoidable source for tensor modes on the brane. Thus it is 
necessary to use the 5D metric-based formalism. This is the subject of the next Section [7] 
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7 Gravitational Wave Perturbations in Brane-World Cos- 
mology 


7.1 Analytical approaches 
The tensor perturbations are given by Equation (273), i.e., (for a flat background brane), 
(ds? = —N?(t, y)dt? + A? (t, y) [bij + fij] da*da? + dy?. (354) 


The transverse traceless fij satisfies Equation (293), which implies, on splitting fij into Fourier 
modes with amplitude f(t, y), 


x À N\. 
i (09-2) 


By the transverse traceless part of Equation (279), the boundary condition is 


1 k? p A N'N, 
na taff b+) (355) 


4 T Vrané = Tij, (356) 
where Tij is the tensor part of the anisotropic stress of matter-radiation on the brane. 

The wave equation (855) cannot be solved analytically except if the background metric functions 
are separable, and this only happens for maximally symmetric branes, i.e., branes with constant 
Hubble rate Ho. This includes the RS case Hp = 0 already treated in Section[2] The cosmologically 
relevant case is the de Sitter brane, Hp > 0. We can calculate the spectrum of gravitational 
waves generated during brane inflation [233], if we approximate slow-roll inflation 
by a succession of de Sitter phases. The metric for a de Sitter brane dS4 in AdS; is given by 


Equations (186, [87] 088) with 


N(t,y) = ny), (357) 
A(t, y) = a(t)n(y), (358) 
n(y) = cosh wy — (1 + ) sinh ply, (359) 
a(t) = ao exp Ho(t — to), (360) 
Hy = po (1 + 2) i (361) 


where u = 47t. 
The linearized wave equation is separable. As before, we separate the amplitude as 
f => em(t)fm(y) where m is the 4D mass, and this leads to: 


. : k? 
Öm + 3Hopm + rn? + “| Pm = 0, (362) 
n! m? 
fin +4— fm + fm = 0. (363) 
n n 
The general solutions for m > 0 are 
3 k 
Pm(t) = exp (— $10") B, (fe) ; (364) 
3/2 p’ 
fin(y) = ry) E372 | 4/1 + Hey (365) 
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Figure 13: Graviton “volcano” potential around the dS4 brane, showing the mass gap. (Figure 
taken from [272].) 


where B, is a linear combination of Bessel functions, Ly /2 is a linear combination of associated 


Legendre functions, and 
fm? 9 


It is more useful to reformulate Equation (363) as a Schrédinger-type equation, 


—V(z)Um = -m Vn, (367) 


using the conformal coordinate 
” dy 1 —ı / Ho 
Z = Zb +f —, zb = — sinh (=) ; (368) 
o ny) Ho H 
and defining Ùm = n3/? fm. The potential is given by (see Figure [[3) 


15H 9 po 
= — 2 + -Hé -3u (1+ —) 6(z- 2), 369 
Ira | 4° u( 7 ) (z — 2p), (369) 


V(2) 


where the last term incorporates the boundary condition at the brane. The “volcano” shape of 
the potential shows how the 5D graviton is localized at the brane at low energies. (Note that 
localization fails for an AdS4 brane [226] [399].) 


The non-zero value of the Hubble parameter implies the existence of a mass gap [155], 


Am = Ho, (370) 


between the zero mode and the continuum of massive KK modes. This result has been generalized: 
For dS4 brane(s) with bulk scalar field, a universal lower bound on the mass gap of the KK tower 


is [149 
[3 
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The massive modes decay during inflation, according to Equation (364), leaving only the zero mode, 
which is effectively a 4D gravitational wave. The zero mode, satisfying the boundary condition 


fo(e,0) = 0, (372) 
is given by 
fo = VEF (Ho/p), (373) 
where the normalization condition so 
2 | \W2\dz = 1 (374) 
Zb 


implies that the function F is given by [274] 


1 1 
=4+,/14+3 
x £ 


-1/2 
F(a) =| VTF- 21n i (375) 


At low energies (Hy < p) we recover the general relativity amplitude: F — 1. At high energies, 
the amplitude is considerably enhanced: 


3H, 
Ho>u > PA Si (376) 
u 


The factor F determines the modification of the gravitational wave amplitude relative to the 


standard 4D result: 
8 (Ho\* 
M3 2r 


The modifying factor F can also be interpreted as a change in the effective Planck mass [149]. 

This enhanced zero mode produced by brane inflation remains frozen outside the Hubble ra- 
dius, as in general relativity, but when it re-enters the Hubble radius during radiation or matter 
domination, it will no longer be separated from the massive modes, since H will not be constant. 
Instead, massive modes will be excited during re-entry. In other words, energy will be lost from the 
zero mode as 5D gravitons are emitted into the bulk, i.e., as massive modes are produced on the 
brane. A phenomenological model of the damping of the zero mode due to 5D graviton emission 
is given in [283]. Self-consistent low-energy approximations to compute this effect are developed 
in [202] [37]. 

At zero order, the low-energy approximation is based on the following [27]. In the 
radiation era, at low energy, the background metric functions obey 


A? = F?(Ho/n). (377) 


A(t, y) > a(t)e ™, N(t,y) > e”. (378) 


To lowest order, the wave equation therefore separates, and the mode functions can be found 
analytically [323] [325] 27]. The massive modes in the bulk, fm(y), are the same as for a Minkowski 
brane. On large scales, or at late times, the mode functions on the brane are given in conformal 
time by 


2 
©) (p) = n72 B j4 ( ZE i 379 
Pm (n) =n a ( a'j (379) 


where ap marks the start of the low-energy regime (pp = A), and B, denotes a linear combination 
of Bessel functions. The massive modes decay on super-Hubble scales, unlike the zero-mode. 
Expanding the wave equation in po/A, one arrives at the first order, where mode-mixing arises. 
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The massive modes yh) (7) on sub-Hubble scales are sourced by the initial zero mode that is 
re-entering the Hubble radius [137]: 


a Ona (1), p2,,(1) 2,3 (1) po 2, (0) 
oe — | em + k*apy + mapi = -47 [mok apò), (380) 


where Imo is a transfer matrix coefficient. The numerical integration of the equations confirms 
the effect of massive mode generation and consequent damping of the zero-mode. 


7.2 Full numerical solutions 


Full numerical solutions for the tensor perturbations have been obtained by the two methods — the 
pseudo-spectral (PS) method used in and the characteristic integration (CI) algorithm 
developed in [70]. It was shown that both methods give identical results and the behaviour of 
gravitational waves on a brane is quite insensitive to the initial conditions in the bulk as for the 
scalar perturbations [381]. Here we summarize the results obtained in [200]. It is convenient to use 
the static bulk metric and consider a moving brane. The simplest initial condition is f(z, T) = const 
when the mode is outside the hubble horizon on the brane. If the brane is static, this would give a 
zero-mode solution f = cos(k(7—79)) where 7 is an initial time. However, due to the motion of the 
brane, which causes the expansion of the brane universe, KK modes are excited and this solution is 
modified. Figure [4] demonstrate this effect. Once the perturbation enters the horizon, non-trivial 
waves are excited in the bulk and the amplitude of the tensor perturbation is damped. Figure [15] 
shows the behaviour of gravitational waves for two different wave numbers. Here ¢, = p/X at the 
time when the mode re-enters the horizon. As for scalar perturbations, we can define the effective 
AD solutions by ignoring the bulk as a reference 


2 
a T (381) 
a(t)? 

href Only takes into account the effect of the high-energy modification of the Friedmann equation. 
Figure [15] shows that the full solution has an additional suppression of the amplitude compared 
with Aref. This suppression is caused by the excitations of KK modes at the horizon crossing as is 
seen in Figure[i4] The suppression is stronger for modes that enter the horizon earlier e, > 1 and 
it becomes negligible at low energies € < 1. 

The ratio |h5p/hrep| evaluated at the low-energy regime long after the horizon re-entry time 
monotonically decreases with the frequency and the suppression of amplitude hsp becomes signif- 
icant above the critical frequency ferit given by 


= 1 crit Aeq 
crit — 5_7 


27E aeg Go 


we e a 
= 0.1 mm 72 km/s - Mpc 3200 l 


This corresponds to a frequency of the mode that enters the horizon when H, = £~' (cf. [203]). 
The ratio |hsp/ħret| obtained from numerical solutions is fitted as 


a (4 -) (383) 


(382) 


hsp 
ret 


with a = 0.76 + 0.01 and 6 = 0.67 + 0.01. 
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Figure 14: The evolution of gravitational waves. We set the comoving wave number to k = v3/£ 


(ex = 1.0). The right panel depicts the projection of the three-dimensional waves of the left panel. 
Figure taken from [200]. 
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Figure 15: Squared amplitude of gravitational waves on the brane in the low-energy (left) and 
the high-energy (right) regimes. In both panels, solid lines represent the numerical solutions. The 


dashed lines are the amplitudes of reference gravitational waves href obtained from Equation (381). 
Figure taken from [200]. 
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There are two important effects on the spectrum in the high-energy regime. Let us first consider 
the non-standard cosmological expansion due to the p?-term. The spectrum of the stochastic 
gravitational waves is modified to 


ee (384) 
fer for f > Terit; 


6w—2 
S pa for f < ferito 
ref — 


where w is an equation of state. This is because gravitational waves re-enter the horizon when the 
p°-term dominates at high frequencies f > ferite In the high-energy radiation dominated phase, 
the spectrum of the stochastic gravitational waves is modified to 


ref xX fee (ferit < fae (385) 


The other effect is the KK-mode excitations. Taking account of the KK-mode excitations, the 


spectrum is calculated as 
2 


h 
EY ee (386) 


href 


w=] 


where we used the fact Qaw œ h? f?. Combining it with the result (883), the spectrum becomes 
nearly flat above the critical frequency: 


Naw « f°, (387) 


which is shown in filled squares in Figure [6] In this figure, the spectrum calculated from the 
reference gravitational waves Qyer is also shown in filled circles. Note that the normalization factor 
of the spectrum is determined as Qgw = 10714 from the CMB constraint. The short-dashed line 
and the solid line represent asymptotic behaviors in the high-frequency region. The spectrum 
taking account of the two high-energy effects seems almost indistinguishable from the standard 
four-dimensional prediction shown in long-dashed line in the figure. In other words, while the 
effect due to the non-standard cosmological expansion enhances the spectrum, the KK-mode effect 
reduces the GW amplitude, which results in the same spectrum as the one predicted in the four- 
dimensional theory. Note that the amplitude taking account of the two effects near f ~% ferit is 
slightly suppressed, which agrees with the results in the previous study for eẹ < 0.3 using the 
Gaussian-normal coordinates [202] discussed in the previous subsection. 

This cancelation of two high energy effects is valid only for w = 1/3. For other equations of 
state, the final spectrum at high frequencies are different from 4D predictions. For example for 
w = 1, New x f? for f > fe while the 4D theory predicts New œx f!. 
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Figure 16: The energy spectrum of the stochastic background of gravitational wave around the 
critical frequency fe in radiation dominated epoch. The filled circles represent the spectrum caused 
by the non-standard cosmological expansion of the universe. Taking account of the KK-mode 
excitations, the spectrum becomes the one plotted by filled squares. In the asymptotic region 
depicted in the solid line, the frequency dependence becomes almost the same as the one predicted 
in the four-dimensional theory (long-dashed line). Figure taken from [200]. 


8 CMB Anisotropies in Brane-World Cosmology 


For the CMB anisotropies, one needs to consider a multi-component source. Linearizing the general 
nonlinear expressions for the total effective energy-momentum tensor, we obtain 


(c= (1 į x $ e=) , 388) 
Po = p+ E (2p+p) +5, 389) 
Got = au (1+ $) +a 390) 
m = Tie (1 oe +) + Ts 391) 
where 
(392) 


p= > Pt p=) Po w= dq 


are the total matter-radiation density, pressure, and momentum density, respectively, and 7, is 
the photon anisotropic stress (neglecting that of neutrinos, baryons, and CDM). 

The perturbation equations in the previous Section[/]form the basis for an analysis of scalar and 
tensor CMB anisotropies in the brane-world. The full system of equations on the brane, including 
the Boltzmann equation for photons, has been given for scalar and tensor |283| perturbations. 
But the systems are not closed, as discussed above, because of the presence of the KK anisotropic 
stress an which acts a source term. 
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In the tight-coupling radiation era, the scalar perturbation equations may be decoupled to give 
an equation for the gravitational potential ®, defined by the electric part of the brane Weyl tensor 
(not to be confused with Env): 


Eu = ViVv®. (393) 
In general relativity, the equation in ® has no source term, but in the brane-world there is a source 


term made up of rË, and its time-derivatives. At low energies (p < åA), and for a flat background 
(K = 0), the equation is [284 


t: 1 2 3 1 
320" + 120! + 2b, = = ne" ang" 4 (= na =) nil (394) 
where x = k/(aH), a prime denotes d/dax, and ©; and që are the Fourier modes of ® and x‘, 
respectively. In general relativity the right hand side is zero, so that the equation may be solved 
for ®,, and then for the remaining perturbative variables, which gives the basis for initializing 


CMB numerical integrations. At high energies, earlier in the radiation era, the decoupled equation 
is fourth order [284]: 


729220 + 38881P + (1782 + 54x2) + 14406), + (90 + x?) h, 


E mi 1 E 1 1 1 2 2 E H 
NF [21 (2) — 0 (ae) ane) (a8 


p x \p x p 
162 2 EN” 4 162 2 E 

_ 30( a 42 + 30(16 “#1 (#)) 398) 
x p x p 


The formalism and machinery are ready to compute the temperature and polarization anisotro- 
pies in brane-world cosmology, once a solution, or at least an approximation, is given for A The 
resulting power spectra will reveal the nature of the brane-world imprint on CMB anisotropies, and 
would in principle provide a means of constraining or possibly falsifying the brane-world models. 
Once this is achieved, the implications for the fundamental underlying theory, i.e., M theory, would 
need to be explored. 

However, the first step required is the solution for Tiy This solution will be of the form 
given in Equation (249). Once G and Fk are determined or estimated, the numerical integration 
in Equation can in principle be incorporated into a modified version of a CMB numerical 
code. The full solution in this form represents a formidable problem, and one is led to look for 
approximations. 


8.1 The low-energy approximation 
The basic idea of the low-energy approximation is to use a gradient 


expansion to exploit the fact that, during most of the history of the universe, the curvature scale 
on the observable brane is much greater than the curvature scale of the bulk (¢ < 1 mm): 


L~ |Ruvop|—/? > £~ |ORagen|- "2 > |V ~ E71 « |dy| ~ £71. (396) 
These conditions are equivalent to the low energy regime, since £? œ A~! and |Ryvag| ~ |Tv: 


(2 
ne A <1. (397) 


Using Equation (896) to neglect appropriate gradient terms in an expansion in (?/L?, the low- 
energy equations can be solved. However, two boundary conditions are needed to determine all 
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functions of integration. This is achieved by introducing a second brane, as in the RS 2-brane 
scenario. This brane is to be thought of either as a regulator brane, whose backreaction on the 
observable brane is neglected (which will only be true for a limited time), or as a shadow brane 
with physical fields, which have a gravitational effect on the observable brane. 

The background is given by low-energy FRW branes with tensions +A, proper times t+, scale 
factors a+, energy densities p+ and pressures p+, and dark radiation densities pg +. The physical 
distance between the branes is ¢d(t), and 


d ed 4 s E J 
i 7 saet H_ =e (H, -å), pe- = elpe +. (398) 


Then the background dynamics is given by 


2 K? 
Hy = += (P4 pe +), (399) 


à e . „2 $ 
d+3H,d—-@ = = [o+ 3p4 + e4(p_ —3p_)]. (400) 


(see for the general background, including the high-energy regime). The dark energy 
obeys pe + = C/ at, where C is a constant. From now on, we drop the +-subscripts which refer to 
the physical, observed quantities. 

The perturbed metric on the observable (positive tension) brane is described, in longitudinal 
gauge, by the metric perturbations ~ and R, and the perturbed radion is d = d+ N. The 
approximation for the KK (Weyl) energy-momentum tensor on the observable brane is 


2 K 


1 
i me (TH +e” T" |) — VEY Ld + UV? d — (veavia+ zava) , (401) 


and the field equations on the observable brane can be written in scalar-tensor form as 


2 a2 Ts 2 1 1 
Gr, = Erpa SOM m 4 L (yey, sev?) + (vx - 5x?) (402) 


3 
x=1-e, wg) =. (403) 


The perturbation equations can then be derived as generalizations of the standard equations. 
For example, the 6G°g equation is [245] 
2d 


© 
i pa E 


1 k2 1 2d = 
Rae — (òp + e=1iöp_) + 


ee 
HY-UR- eG aaa 


2 
e-* — 1 3 
= -— [(1- H) + (i- H)? wus onde — ik- oN l 


(404) 


The trace part of the perturbed field equation shows that the radion perturbation determines the 
crucial quantity, d7¢: 


2 
R +4 = —-——N = -ka ne, (405) 
e2d — | 
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where the last equality follows from Equation (825). The radion perturbation itself satisfies the 
wave equation 


+ Ch A . x aa k2 
Ñ + (3H — 2d) Ñ — (2H +4H? + 2} - 6Hå— 2d) N + SN 
a 


hb + 3dR 4 ( 2d 6Hd + 2d”) w 


2 


K 
6 
P4 


[50 — 36p + e-24(5p_ — 36p-)| (406) 


A new set of variables 


.,# turns out be very useful 247]: 


| 
2%) - Hee Ee, 


k2 3 
a2 oe . 
= -9+ — yz (Ë + 2HE), 


g~ a 
N =- p4} pate. (407) 
Equation (405) gives 
z ad \ . 1k? Qe2d k? od 


The variable Æ determines the metric shear in the bulk, whereas y+ give the brane displace- 
ments in transverse traceless gauge. The latter variables have a simple relation to the curvature 
perturbations on large scales (restoring the +-subscripts): 


H? (ọọ 
arscgpr E fE o I, 409 
Cot p+ + A (= +s) (409) 
where f+ = df+/dt+. 
8.2 The simplest model 
The simplest model is the one in which 
pe=0=d (410) 


in the background, with p_/p_ = p/p. The regulator brane is assumed to be far enough away that 
its effects on the physical brane can be neglected over the timescales of interest. By Equation (400) 
it follows that 

p- = —pe4 (411) 


i.e., the matter on the regulator brane must have fine-tuned and negative energy density to prevent 
the regulator brane from moving in the background. With these assumptions, and further assuming 
adiabatic perturbations for the matter, there is only one independent brane-world parameter, i.e., 
the parameter measuring dark radiation fluctuations: 


5 
i=, (412) 
Prad 


This assumption has a remarkable consequence on large scales: The Weyl anisotropic stress 
ne terms in the Sachs—Wolfe formula (827) cancel the entropy perturbation from dark radiation 
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Figure 17: The CMB power spectrum with brane-world effects, encoded in the dark radiation 
fluctuation parameter dC, as a proportion of the large-scale curvature perturbation for matter 


(denoted Ç» in the plot). (Figure taken from [247].) 


fluctuations, so that there is no difference on the largest scales from the standard general relativity 
power spectrum. On small scales, beyond the first acoustic peak, the brane-world corrections are 
negligible. On scales up to the first acoustic peak, brane-world effects can be significant, changing 
the height and the location of the first peak. These features are apparent in Figure[17] However, it 
is not clear to what extent these features are general brane-world features (within the low-energy 
approximation), and to what extent they are consequences of the simple assumptions imposed on 


the background. Further work remains to be done. 
A related low-energy approximation, using the moduli space approximation, has been developed 


for certain 2-brane models with bulk scalar field [364] BI]. The effective gravitational action on 


the physical brane, in the Einstein frame, is 


E 1 4 12a? 2 6 2 


where a is a coupling constant, and ¢ and y are moduli fields (determined by the zero-mode of 
the bulk scalar field and the radion). Figure [I8] shows how the CMB anisotropies are affected by 


the x-field. 
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Figure 18: The CMB power spectrum with brane-world moduli effects from the field x in Equa- 
tion (413). The curves are labelled with the initial value of x. (Figure taken from BI.) 
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9 DGP Models: Modifying Gravity at Low Energies 
9.1 ‘Self-accelerating’ DGP 


Most brane-world models modify general relativity at high energies. The Randall-Sundrum models 
discussed up to now are a typical example. At low energies, H£ « 1, the zero-mode of the graviton 
dominates on the brane, and general relativity is recovered to a good approximation. At high 
energies, Hf >> 1, the massive modes of the graviton dominate over the zero-mode, and gravity 
on the brane behaves increasingly 5-dimensional. On the unperturbed FRW brane, the standard 
energy-conservation equation holds, but the Friedmann equation is modified by an ultraviolet 
correction, (Gép)*. At high energies, gravity “leaks” off the brane and H? « p?. 

By contrast, the brane-world model of Dvali-Gabadadze—Porrati (DGP), which was gen- 
eralized to cosmology by Deffayet [116], modifies general relativity at low energies. This model 
produces ‘self-acceleration’ of the late-time universe due to a weakening of gravity at low energies. 
Like the RS model, the DGP model is a 5D model with infinite extra dimensions? 

The action is given by 


= B da y —g®) RO) +f 
bulk 


167G Te brane 


dtr yg r| l (414) 


The bulk is assumed to be 5D Minkowski spacetime. Unlike the AdS bulk of the RS model, the 
Minkowski bulk has infinite volume. Consequently, there is no normalizable zero-mode of the 4D 
graviton in the DGP brane-world. Gravity leaks off the 4D brane into the bulk at large scales, 
r © re, where the first term in the sum dominates. On small scales, gravity is effectively 
bound to the brane and 4D dynamics is recovered to a good approximation, as the second term 
dominates. The transition from 4D to 5D behaviour is governed by the crossover scale re. For a 
Minkowski brane, the weak-field gravitational potential behaves as 


rol forr & re 
á fra forr > re ` (415) 


On a Friedmann brane, gravity leakage at late times in the cosmological evolution can initiate 
acceleration — not due to any negative pressure field, but due to the weakening of gravity on the 
brane. 4D gravity is recovered at high energy via the lightest massive modes of the 5D graviton, 
effectively via an ultra-light meta-stable graviton. 

The energy conservation equation remains the same as in general relativity, but the Friedmann 
equation is modified [116]: 


p+3H(p+p)=0, (416) 
K 1 K 8rG 


The Friedmann equation can be derived from the junction condition or the Gauss-Codazzi equation 
as in the RS model. To arrive at Equation (417) we have to take a square root which implies a 
choice of sign. As we shall see, the above choice has the advantage of leading to acceleration — 
but at the expense of a ‘ghost’ (negative kinetic energy) mode in the scalar graviton sector. The 
‘normal’ (non-self-accelerating) DGP model, where the opposite sign of the square root is chosen, 
has no ghost, and is discussed below. 

From Equation we infer that at early times, i.e., Hre >> 1, the general relativistic 
Friedman equation is recovered. By contrast, at late times in an expanding CDM universe, with 


3DGP brane-worlds without Z2 symmetry have also been considered; see e.g. [396]. 
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Figure 19: Joint constraints [solid thick (blue)] from the SNLS data [solid thin (yellow)], the 
BAO peak at z = 0.35 [dotted (green)] and the CMB shift parameter from WMAP3 [dot-dashed 
(red)]. The left plot shows LCDM models, the right plot shows DGP. The thick dashed (black) 
line represents the flat models, Qg = 0. (From [308].) 
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Figure 20: The constraints from SNe and CMB/BAO on the parameters in the DGP model. The 
flat DGP model is indicated by the vertical dashed-dotted line; for the MLCS light curve fit, the 
flat model matches to the date very well. The SALT-II light curve fit to the SNe is again shown 
by the dotted contours. The combined constraints using the SALT-II SNe outlined by the dashed 
contours represent a poorer match to the CMB/BAO for the flat model. (From [404].) 
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p «xa? — 0, we have 
1 
H = H% = —, (418) 


Te 


so that expansion accelerates and is asymptotically de Sitter. The above equations imply 


1 


_ K 
B= = —4rGp |1 + Se (419) 
a y 1 + 32rGr2p/3 
In order to achieve self-acceleration at late times, we require 
re 2 Hg', (420) 


since Ho < Ho. This is confirmed by fitting supernova observations, as shown in Figures [9] and 
The dimensionless cross-over parameter is defined as 


1 
Qr. = ——— 421 
c 4( Hore)? J ( ) 
and the LCDM relation, 
Am +a +Nk =1, (422) 


is modified to 
Om + 2fOr,f1-QK+F0K =1. (423) 


LCDM and DGP can both account for the supernova observations, with the fine-tuned values 
A ~ HÊ and re ~ Ho 1 respectively. When we add further constraints to the expansion history 
from the baryon acoustic oscillation peak at z = 0.35 and the CMB shift parameter, the DGP flat 
models are in strong tension with the data, whereas LCDM models provide a consistent fit. This 
is evident in Figures [19] and 20] though this conclusion depends on a choice of light curve fitters 
in the SNe observations. The open DGP models provide a somewhat better fit to the geometric 
data — essentially because the lower value of Qm favoured by supernovae reduces the distance to 
last scattering and an open geometry is able to extend that distance. For a combination of SNe, 
CMB shift and Hubble Key Project data, the best-fit open DGP also performs better than the flat 
DGP [405], as shown in Figure 

Observations based on structure formation provide further evidence of the difference between 
DGP and LCDM, since the two models suppress the growth of density perturbations in different 
ways [297]. The distance-based observations draw only upon the background 4D Friedman 
equation in DGP models — and therefore there are quintessence models in general relativity 
that can produce precisely the same supernova distances as DGP. By contrast, structure formation 
observations require the 5D perturbations in DGP, and one cannot find equivalent quintessence 
models [253]. One can find 4D general relativity models, with dark energy that has anisotropic 
stress and variable sound speed, which can in principle mimic DGP [265]. However, these models 
are highly unphysical and can probably be discounted on grounds of theoretical consistency. 

For LCDM, the analysis of density perturbations is well understood. For DGP the perturbations 
are much more subtle and complicated [253]. Although matter is confined to the 4D brane, gravity 
is fundamentally 5D, and the 5D bulk gravitational field responds to and back-reacts on 4D density 
perturbations. The evolution of density perturbations requires an analysis based on the 5D nature 
of gravity. In particular, the 5D gravitational field produces an effective “dark” anisotropic stress 
on the 4D universe, as discussed in Section 3.4] If one neglects this stress and other 5D effects, 
and simply treats the perturbations as 4D perturbations with a modified background Hubble rate 
— then as a consequence, the 4D Bianchi identity on the brane is violated, i.e., V’Gy, # 0, and the 
results are inconsistent. When the 5D effects are incorporated [I], the 4D Bianchi identity 
is automatically satisfied. (See Figure [22]) 

There are three regimes governing structure formation in DGP models: 
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Figure 21: The difference in x? between best-fit DGP (flat and open) and best-fit (flat) LCDM, 
using SNe, CMB shift and Hp Key Project data. (From [405].) 
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On small scales, below the Vainshtein radius (which for cosmological purposes is roughly the 
scale of clusters), the spin-0 scalar degree of freedom becomes strongly coupled, so that the 
general relativistic limit is recovered [258]. 


On scales relevant for structure formation, i.e., between cluster scales and the Hubble ra- 
dius, the spin-0 scalar degree of freedom produces a scalar-tensor behaviour. A quasi-static 
approximation (as in the Newtonian approximation in standard 4D cosmology) to the 5D 
perturbations shows that DGP gravity is like a Brans—Dicke theory with parameter [253 


wap = 5(8 I), (424) 


a ie H 2K 
2 2 


At late times in an expanding universe, when Hr, Z 1, it follows that 6 < 1, so that wep < 0. 


~N 


(This is a signal of the ghost pathology in DGP, which is discussed below.) 


Although the quasi-static approximation allows us to analytically solve the 5D wave equation 
for the bulk degree of freedom, this approximation breaks down near and beyond the Hubble 
radius. On super-horizon scales, 5D gravity effects are dominant, and we need to solve 
numerically the partial differential equation governing the 5D bulk variable [7I]. 


Figure 22: The growth factor g(a) = A(a)/a for LCDM (long dashed) and DGP (solid, thick), as 
well as for a dark energy model with the same expansion history as DGP (short dashed). DGP-4D 
(solid, thin) shows the incorrect result in which the 5D effects are set to zero. (From [253].) 


On subhorizon scales relevant for linear structure formation, 5D effects produce a difference 
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Figure 23: Left: Numerical solutions for DGP density and metric perturbations, showing also 
the quasi-static solution, which is an increasingly poor approximation as the scale is increased. 
(From [71].) Right: Constraints on DGP (the open model in Figure2I]that provides a best fit to 
geometric data) from CMB anisotropies (WMAP5). The DGP model is the solid curve, QCDM 
(short-dashed curve) is the GR quintessence model with the same background expansion history 
as the DGP model, and LCDM is the dashed curve (a slightly closed model that gives the best fit 
to WMAP5, HST and SNLS data). (From [143].) 
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between ¢ and —w [253]: 


1 
k?¢ = 40Ga? (1 - 5) pA, (426) 
1 
kw = —41Ga? (1 + 5) pA, (427) 
so that there is an effective dark anisotropic stress on the brane: 
8rGa? 
kele +) = 5 pA. (428) 
36 
The density perturbations evolve as 
y : 1 
A+2HA -4G (1-73) =o. (429) 


The linear growth factor, g(a) = A(a)/a (i.e., normalized to the flat CDM case, A x a), is shown 
in Figure This illustrates the dramatic suppression of growth in DGP relative to LCDM — 
from both the background expansion and the metric perturbations. If we parametrize the growth 
factor in the usual way, we can quantify the deviation from general relativity with smooth dark 


energy [295]: 


0.55 + 0.05[1 + w(z = 1)] GR, smooth DE 
= = y a 
Ta A. T a DGP oo 


Observational data on the growth factor [189| are not yet precise enough to provide meaningful 
constraints on the DGP model. Instead, we can look at the large-angle anisotropies of the CMB, 
i.e., the ISW effect. This requires a treatment of perturbations near and beyond the horizon 
scale. The full numerical solution has been given by [71], and is illustrated in Figure 23] The 
CMB anisotropies are also shown in Figure|23} as computed in using a scaling approximation 
to the super-Hubble modes [379] (the accuracy of the scaling ansatz is verified by the numerical 
results [71]). 

It is evident from Figure 23] that the DGP model, which provides a best fit to the geometric 
data (see Figure BIJ, is in serious tension with the WMAP5 data on large scales. The problem 
arises because there is a large deviation of ¢_- = (@ — w)/2 in the DGP model from the LCDM 
model. This deviation, i.e., a stronger decay of d_, leads to an over-strong ISW effect (which is 
determined by b_), in tension with WMAP5 observations. 

As a result of the combined observations of background expansion history and large-angle 
CMB anisotropies, the DGP model provides a worse fit to the data than LCDM at about the 50 
level [143]. Effectively, the DGP model is ruled out by observations in comparison with the LCDM 
model. 

In addition to the severe problems posed by cosmological observations, a problem of theo- 
retical consistency arises from the fact that the late-time asymptotic de Sitter solution in DGP 
cosmological models has a ghost. The ghost is signaled by the negative Brans—Dicke parameter 
in the effective theory that approximates the DGP on cosmological subhorizon scales: The ex- 
istence of the ghost is confirmed by detailed analysis of the 5D perturbations in the de Sitter 
limit [248]. The DGP ghost is a ghost mode in the scalar sector of the gravitational 
field — which is more serious than the ghost in a phantom scalar field. It effectively rules out the 
DGP, since it is hard to see how an ultraviolet completion of the DGP can cure the infrared ghost 
problem. However, the DGP remains a valuable toy model for illustrating the kinds of behaviour 
that can occur from a modification to Einstein’s equations — and for developing cosmological tools 
to test modified gravity and Einstein’s theory itself. 
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9.2 ‘Normal’ DGP 


The self-accelerating DGP is effectively ruled out as a viable cosmological model by observations 
and by the problem of the ghost in the gravitational sector. Indeed, it may be the case that 
self-acceleration generically comes with the price of ghost states. The ‘normal’ (i.e., non-self- 
accelerating and ghost-free) branch of the DGP [374], arises from a different embedding of the 
DGP brane in the Minkowski bulk (see Figure 24). In the background dynamics, this amounts 
to a replacement re —> —re in Equation (417) — and there is no longer late-time self-acceleration. 
Therefore, it is necessary to include a A term in order to accelerate the late universe: 


K 1] K  8rnG A 
a 24 = = —_— =. 431 
salar ae oe 3 PTR gal) 


(Normal DGP models with a quintessence field have also been investigated [88].) Using the di- 
mensionless crossover parameter defined in Equation (421), the densities are related at the present 


time by 
VJ1—-2 = —V Qr, + Qr, + Xm + Da, (432) 


which can be compared with the self-accelerating DGP relation (423). 
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Figure 24: Left: The embedding of the self-accelerating and normal branches of the DGP brane 
in a Minkowski bulk. (From [81].) Right: Joint constraints on normal DGP (flat, K = 0) from 
SNLS, CMB shift (WMAP3) and BAO (z = 0.35) data. The best-fit is the solid point, and is 
indistinguishable from the LCDM limit. The shaded region is unphysical and its upper boundary 
represents flat LCDM models. (From [281].) 


An interesting feature of the normal branch is the ‘degravitation’ property, i.e., that A is 
effectively screened by 5D gravity effects. This follows from rewriting the modified Friedmann 
equation (431) in standard general relativistic form, with 


K 
brste ee Sh, (433) 
Te a 


Thus, 5D gravity in normal DGP can in principle reduce the bare vacuum energy significantly [374] 
[300]. However, Figure[24]shows that best-fit flat models, using geometric data, only admit insignif- 
icant screening [281]. The closed models provide a better fit to the data [168], and can allow a 
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bare vacuum energy term with Qa > 1, as shown in Figure 5) This does not address the funda- 
mental problem of the smallness of Q4, but is nevertheless an interesting feature. We can define 
an effective equation-of-state parameter via 


Aeg + 3H(1 + Wert ) Neg = 06 (434) 
At the present time (setting K = 0 for simplicity), 


(Q +A ) <1 


Mamta ~~” (435) 


Weff,0 = =] = 


where the inequality holds because Qm < 1. This reveals another important property of the nor- 
mal DGP model: effective phantom behaviour of the recent expansion history. This is achieved 
without any pathological phantom field (similar to what can be done in scalar-tensor theories ). 
Furthermore, there is no “big rip” singularity in the future associated with this phantom acceler- 
ation, unlike the situation that typically arises with phantom fields. The phantom behaviour in 
the normal DGP model is also not associated with any ghost problem — indeed, the normal DGP 
branch is free of the ghost that plagues the self-accelerating DGP [81I]. 


a 1.054 
Wes 
> 1}@ ETL 
< 0.954 
wy 
g 0.94 — quasi-static approximation 
S ogi ° k= 0.0005 h Mpc! 
a e k = 0.002 h Mpc™' DN 
Z 084 o k=0.01hMpc™ KAN 
So 0.754 
0 ca RTT BESS Ee Ke Be EE Fee Rae a 
—~ -0.02 4 
> 
26:00 + 0.044 
= -0.06 4 t go 
7 _0.08 | — quasi-static approximation o i 
$ + k = 0.0005 h Mpc”! S 
+ nDGP2 eo 1] 2 k=0.002h Mpc! 4 
= 0.12] = k=0.01hMpc™* 
0.05 e nDGP3 : 
Š 1.054 
1+ ee - ---- 
0.10 | 0.954 
i Zia 0.94 — quasi-static approximation 2 
0854 -+ k= 0.0005% Mpc”! o 
| e k = 0.002 h Mpc™' iy 
t 1 @ 08] a k=0.01AMpc7! [i 
-0.15 1 La 1 La . L 1 0.754 e=—l4 
9:3 ne 13 a0 2:5 2-18 -16 -1.4 -1.2 -1 -0.8 -0.6 -0.4 02 0 


scale factor loga 


Figure 25: Left: Joint constraints on normal DGP from SNe Gold, CMB shift (WMAP3) and Ho 
data in the projected curvature-A plane, after marginalizing over other parameters. The best-fits 
are the solid points, corresponding to different values of Qm. (From [168].) Right: Numerical 
solutions for the normal DGP density and metric perturbations, showing also the quasistatic so- 
lution, which is an increasingly poor approximation as the scale is increased. Compare with the 
self-accelerating DGP case in Figure[23] (From [71].) 
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Perturbations in the normal branch have the same structure as those in the self-accelerating 
branch, with the same regimes — i.e., below the Vainshtein radius (recovering a GR limit), up to the 
Hubble radius (Brans—Dicke behaviour), and beyond the Hubble radius (strongly 5D behaviour). 
The quasistatic approximation and the numerical integrations can be simply repeated with the 
replacement re — —re (and the addition of A to the background). In the sub-Hubble regime, 
the effective Brans—Dicke parameter is still given by Equations and (425), but now we have 
wpgp > 0-and this is consistent with the absence of a ghost. Furthermore, a positive Brans—Dicke 
parameter signals an extra positive contribution to structure formation from the scalar degree of 
freedom, so that there is less suppression of structure formation than in LCDM - the reverse of 
what happens in the self-accelerating DGP. This is confirmed by computations, as illustrated in 
Figure 

The closed normal DGP models fit the background expansion data reasonably well, as shown 
in Figure 25) The key remaining question is how well do these models fit the large-angle CMB 
anisotropies, which is yet to be computed at the time of writing. The derivative of the ISW potential 
ġ— can be seen in Figure and it is evident that the ISW contribution is negative relative to 
LCDM at high redshifts, and goes through zero at some redshift before becoming positive. This 
distinctive behaviour may be contrasted with the behaviour in f(R) models (see Figure 26): both 
types of model lead to less suppression of structure than LCDM, but they produce different ISW 
effects. However, in the limit r, — oo, normal DGP tends to ordinary LCDM, hence observations 
which fit LCDM will always just provide a lower limit for re. 
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Figure 26: Top: Measurement of the cross-correlation functions between six different galaxy data 
sets and the CMB, reproduced from [167]. The curves show the theoretical predictions for the 
ISW-galaxy correlations at each redshift for the LCDM model (black, dashed) and the three nDGP 
models, which describe the 1—o region of the geometry test from Figure25] (From [168].) Bottom: 
Theoretical predictions for a family of f (R) theqyses compared with the ISW data measuring 
the angular CCF between the CMB and six galaxy catalogues. The model with Bo = 0 is equivalent 
to LCDM, while increasing departures from GR produce negative cross-correlations. (From [{166].) 


10 6-Dimensional Models 


For brane-world models in 6-dimensional spacetime, the codimension of a brane is two and the 
behaviour of gravity is qualitatively very different from the codimension one brane-world models. 
Here we briefly discuss some important examples and features of 6-dimensional models. 


10.1 Supersymmetric Large Extra Dimensions (SLED) Model 


This is a supersymmetric version of the Einstein-Maxwell model (see for 
reviews). The bosonic part of the action is given by 


1 1 

S= few — (R — du o0™ e) — ge Fun FN — e? Ng ; (436) 
6 

There exists a solution where the dilaton ¢ is constant, 6 = ¢9. The solution for the 6D spacetime 


is given by MRE 
ds? = Nuv dx da” + yigdu' dx). (437) 


The gauge field is taken to consist of magnetic flux threading the extra dimensional space so that 
the field strength takes the form 


Fij = V1Bokiy; (438) 


where Bo is a constant, y is the determinant of yij and €;; is the antisymmetric tensor normalized 
as €12 = 1. All other components of Fy» vanish. A static and stable solution is obtained by choosing 
the extra-dimensional space to be a two-sphere 


yijdx’dx = a(d?’ + sin? Ody’). (439) 
The magnetic field strength Bo and the radius ao are fixed by the cosmological constant 
(440) 
The constant value ġo is determined by the condition that the potential for o has minimum [154 
1 
V'(¢0) = E Boe” + Age® =0. (441) 
This is exactly the condition to have a flat geometry on the brane (see Equation (440)) 
Bre = 2Ace*, (442) 


Thus without tuning the 6D cosmological constant, it is possible to obtain a flat 4D spacetime. 
Now we add branes to this solution [72]. The brane action is given by 


S4 = —0 i diry z7. (443) 
The solution for the extra dimensions is now given by 
yijdzřdzÍ = a2(db? + a? sin? Ody’), (444) 
where a 
a=1- 08 = 3h toe (445) 
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The coordinate y ranges from 0 to 27. Thus the effect of the brane makes a deficit angle 6 = 
2n(1 — a) in the bulk (see Figure 27). This is a 6D realization of the ADD model including the 
self-gravity of branes. An interesting property of this model is that regardless of the tension of 
the brane, the 4D specetime on the brane is flat. Thus this could solve the cosmological constant 
problem - any vacuum energy on the brane only changes a geometry of the extra-dimensions and 
does not curve the 4D spacetime. This idea of solving the cosmological constant problem is known 
as self-tuning. 


Figure 27: Removing a wedge from a sphere and identifying opposite sides to obtain a rugby ball 
geometry. Two equal-tension branes with conical deficit angles are located at either pole; outside 
the branes there is constant spherical curvature. From [72]. 


We should note that there have been several objections to the idea of self-tuning [422] [423]. 
Consider that a phase transition occurs and the tension of the brane changes from gı to o2. 
Accordingly, a changes from a, = 1 — 0;/(27M@) to a2 = 1 — o2/(27M@). The magnetic flux is 
conserved as the gauge field strength is a closed form, dF = 0. Then the magnetic flux which is 
obtained by integrating the field strength over the extra dimensions should be conserved 


Dp = 4ra Bo, = 4ra Bo,2. (446) 


The relation between Ag and Bo, Equation (440), that ensures the existence of Minkowski branes 
cannot be imposed both for Bo = Bo,ı and Bo = Bo,2 when a; Æ a2 unless o changes. Moreover, 
the quantization condition must be imposed on the flux ®g. What happens is that a modulus, 
which is a combination of ¢ and the radion describing the size of extra-dimension, acquires a 
runaway potential and the 4D spacetime becomes non-static. 

An unambiguous way to investigate this problem is to study the dynamical solutions directly in 
the 6D spacetime. However, once we consider the case where the tension becomes time dependent, 
we encounter a difficulty to deal with the branes |422|. This is because for co-dimension 2 branes, 
we encounter a divergence of metric near the brane if we put matter other than tension on a brane. 
Hence, without specifying how we regularize the branes, we cannot address the question what 
will happen if we change the tension. Is the self-tuning mechanism at work and does it lead to 
another static solution? Or do we get a dynamical solution driven by the runaway behaviour of 
the modulus field? 

There was a negative conclusion on the self-tuning in this supersymmetric model for a particular 
kind of regularization [423]. However, the answer could depend on the regularization of 
branes and the jury remains out. It is important to study the time-dependent dynamics in the 6D 
spacetime and the regularization of the branes in detail [413] [60] [62] AA [61] [33] [111] [341] [353] (100). 
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10.2 Cosmology in 6D brane-world models 


It is much harder to obtain cosmological solutions in 6D models than in 5-dimensional models. In 
5D brane-world models, cosmological solutions can be obtained by considering a moving brane in a 
static bulk spacetime. This is because the motion of the homogeneous and isotropic brane does not 
change the bulk spacetime thanks to Birkoff’s theorem. However, this is no longer the case in 6D 
spacetime. Thus we need to find time dependent bulk solutions that are coupled to a motion of a 
brane, which requires us to solve non-linear partial differential equations numerically. Moreover, as 
is mentioned above, there appears a curvature singularity if one considers an infinitely thin brane 
and puts matter other than tension on the brane [423]. Thus we need some regularization 
scheme to find cosmological solutions. 

One of the popular ways to regularize a brane is to promote a brane that is a point-like object 
in two extra-dimensions to a 5D ring [356]. The ring is a codimension one object and it is possible 
to consider a motion of this ring. However, a problem is that this ring brane is not homogeneous 
as one of dimensions on the brane is compact and this breaks Birkhoff’s theorem. It has been 
shown that cosmology obtained by the motion of a brane in a static bulk shows pathological 
behaviour [354] [328]. This indicates that we need to find fully time-dependent bulk solutions. 

At the moment, the only accessible way is to solve the 6D bulk spacetime using the gradient 
expansion methods assuming physical scales on a brane are much lower than mass scales in the 
bulk [234]. It is still an open question what are high energy effects in 6D models. 

As in 5D models, there are many generalizations such as the inclusion of the induced gravity 
term on a brane and the Gauss—Bonnet term in the bulk [82]. 


10.3 Cascading brane-world model 


This is a 6D extension of the DGP model 112). The action is given by 


1 1 
g= One Pa /—6 G9) R+ One Pr /—8gOR 


4—brane 


1 TH 
+ / d‘x =g [Or ar Lmatter i (447) 
3—brane 


-2 
2K 


As in the 5D DGP model, we can define the cross over scale. In this model, there are two cross 
over scales: 


Mè 
Or, = ve On = ve (448) 
Assuming that r, > Öre, the gravitational potential on the 3-brane cascades from a 1 /r (4D 
gravity) regime at short scales, to a 1/r? (5D gravity) regime at intermediate distance and finally 
a 1/r? (6D gravity) regime at large distances. 

This model addresses several fundamental issues in induced gravity models in 6D spacetime. 
Without the induced gravity term on the 4-brane, the 6D graviton propagator diverges logarith- 
mically near the position of the 3-brane [163]. On the other hand, the graviton propagator on 
the 3-brane in this model behaves like G(p) > log(p® re) in the Ms — 0 limit where p is 4D 
momentum. Thus the cross-over scale r, acts as a cut-off for the propagator and it remains 
finite even at the position of the 3-brane. 

A more serious issue in the induced gravity model is that most constructions seem to be plagued 
by ghost instabilities [129] [153]. Usually, the regularization of the brane is necessary and it depends 
on the regularization scheme whether there appears a ghost or not [244]. In the cascading model, 
there is still a ghost if the 3-brane has no tension. However, it was shown that by adding a tension 
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to the 3-brane, this ghost disappears 113}. More precisely, there is a critical tension 


crit = Oram (449) 
and above the critical tension À > Aecrit, the model is free of ghosts. It is still unclear what is the 
meaning of having the critical tension for the existence of the ghost. For example, what happens 
if there is a phase transition on the 3-brane and the tension changes from A > Acrit to A < Acrit? 
This remains a very interesting open question. This model also provides an interesting insight 
into the “degravitation” mechanism by which the cosmological constant does not gravitate on the 
3-brane [133] [12]. 

Cosmological solutions in the cascading brane model are again notoriously difficult to find be- 
cause it is necessary to find 6-dimensional solutions that depend on time and two extra-coordinates [I]. 
The simplest de Sitter solutions have been obtained [326]. Interestingly, there exists a self- 
accelerating solution for the 3-brane even when the solution for the 4-brane is in the normal 
branch. It is still not clear whether this self-accelerating solution is stable or not and it is crucial 
to check the stability of this new self-accelerating solution. 

A similar class of models includes intersecting branes [103]. In this model, we have 
two 4-branes that intersect and a 3-brane sits at the intersection. Again there are self-accelerating 
de Sitter solutions and cosmology has been studied by considering a motion of one of the 4-branes. 
A model without a 4-brane has been studied by regularizing a 3-brane by promoting it to a 5D 


ring brane [220] [219]. 
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11 Conclusion 


Simple brane-world models provide a rich phenomenology for exploring some of the ideas that 
are emerging from M theory. The higher-dimensional degrees of freedom for the gravitational 
field, and the confinement of standard model fields to the visible brane, lead to a complex but 
fascinating interplay between gravity, particle physics, and geometry, that enlarges and enriches 
general relativity in the direction of a quantum gravity theory. 

This review has attempted to show some of the key features of brane-world gravity from the 
perspective of astrophysics and cosmology, emphasizing a geometric approach to dynamics and 
perturbations. It has focused mainly on 1-brane RS-type brane-worlds, but also considered the 
DGP brane-world models. The RS-type models have some attractive features: 


e They provide a simple 5D phenomenological realization of the Horava—Witten supergravity 
solutions in the limit where the hidden brane is removed to infinity, and the moduli effects 
from the 6 further compact extra dimensions may be neglected. 


e They develop a new geometrical form of dimensional reduction based on a strongly curved 
(rather than flat) extra dimension. 


e They provide a realization to lowest order of the AdS/CFT correspondence. 
e They incorporate the self-gravity of the brane (via the brane tension). 


e They lead to cosmological models whose background dynamics are completely understood 
and reproduce general relativity results with suitable restrictions on parameters. 


The review has highlighted both the successes and some remaining open problems of the RS 
models and their generalizations. The open problems stem from a common basic difficulty, i.e., 
understanding and solving for the gravitational interaction between the bulk and the brane (which 
is nonlocal from the brane viewpoint). The key open problems of relevance to astrophysics and 
cosmology are 


e to find the simplest realistic solution (or approximation to it) for an astrophysical black hole 
on the brane, and settle the questions about its staticity, Hawking radiation, and horizon; 
and 


e to develop realistic approximation schemes (building on recent work [401| [390} [402] [403] 
[247] [364] [51] [202] [137] ) and manageable numerical codes (building on [247 [B64 [51] [202 [137] ) 
to solve for the cosmological perturbations on all scales, to compute the CMB anisotropies 
and large-scale structure, and to impose observational constraints from high-precision data. 


The RS-type models are the simplest brane-worlds with curved extra dimension that allow 
for a meaningful approach to astrophysics and cosmology. One also needs to consider gener- 
alizations that attempt to make these models more realistic, or that explore other aspects of 
higher-dimensional gravity which are not probed by these simple models. Two important types of 
generalization are the following: 


e The inclusion of dynamical interaction between the brane(s) and a bulk scalar field, so that 
the action is 
1 K 
S=53 Par/—Og [Or — 62049046 — 2A5(®)| + | dryg ao) + =z + Lmatter 
Ks brane(s) v5 
(450) 


91 


(see 
B27 222, B17) [150] [18]). The scalar field could represent a bulk dilaton of the gravitational 
sector, or a modulus field encoding the dynamical influence on the effective 5D theory of an 
extra dimension other than the large fifth dimension 
180). 


For two-brane models, the brane separation introduces a new scalar degree of freedom, the 
radion. For general potentials of the scalar field which provide radion stabilization, 4D 
Einstein gravity is recovered at low energies on either brane [285]. (By contrast, in 
the absence of a bulk scalar, low energy gravity is of Brans—Dicke type [156].) In particular, 
such models will allow some fundamental problems to be addressed: 


— The hierarchy problem of particle physics. 


— An extra-dimensional mechanism for initiating inflation (or the hot radiation era with 
super-Hubble correlations) via brane interaction (building on the initial work in [34] 


[221] [230] [216] [342] [406} 275] [319] [415] R6] 98] [301] [40] [157] [158]). 


— An extra-dimensional explanation for the dark energy (and possibly also dark matter) 
puzzles: Could dark energy or late-time acceleration of the universe be a result of 
gravitational effects on the visible brane of the shadow brane, mediated by the bulk 
scalar field? 


e The addition of stringy and quantum corrections to the Einstein-Hilbert action, including 
the following: 


— Higher-order curvature invariants, which arise in the AdS/CFT correspondence as next- 
to-leading order corrections in the CFT. The Gauss—Bonnet combination in particular 
has unique properties in 5D, giving field equations which are second-order in the bulk 
metric (and linear in the second derivatives), and being ghost-free. The action is 


S= = Pay/—Og [OR — 24s +o OR? 4 ORG ORAP + ORawoy ORABCP) 
5 


+ J dryg |- + = + Liter} (451) 
brane Kg 
where a is the Gauss-Bonnet coupling constant, related to the string scale. The cos- 
mological dynamics of these brane-worlds is investigated in [122] [346] [348] [344] (762) 
[291] [37] [820] [293] [82] [126] [19] 125] [79] [B12]. In [20] it is shown that the black string 
solution of the form of Equation is ruled out by the Gauss—Bonnet term. In this 
sense, the Gauss—Bonnet correction removes an unstable and singular solution. 
In the early universe, the Gauss-Bonnet corrections to the Friedmann equation have 
the dominant form 

H? x p?/3 (452) 


at the highest energies. If the Gauss-Bonnet term is a small correction to the Einstein- 
Hilbert term, as may be expected if it is the first of a series of higher-order corrections, 
then there will be a regime of RS-dominance as the energy drops, when H? « p°. Finally 
at energies well below the brane tension, the general relativity behaviour is recovered. 


— Quantum field theory corrections arising from the coupling between brane matter and 
bulk gravitons, leading to an induced 4D Ricci term in the brane action. The original 
induced gravity brane-world is the DGP model , which we investigated 
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in this review as an alternative to the RS-type models. Another viewpoint is to see the 
induced-gravity term in the action as a correction to the RS action: 


í K 
c= 5a J drey -64g |r = 24s] + 1 diry Zg før —A+ a Loarel (53) 
5 b 5 


rane 


where £ is a positive coupling constant. 


The cosmological models have been analyzed in [118] [239] [27] [231] B75] [394] [373] 
A BI 299 [340] [87] 241]. (Brane-world black holes with induced gravity are 
investigated in [242].) Unlike RS-type models, DGP models lead to 5D behaviour on 
large scales rather than small scales. Then on an FRW brane, the late-universe 5D 
behaviour of gravity can naturally produce a late-time acceleration, even without dark 
energy, although the self-accelerating models suffer from a ghost. Nevertheless, the 
DGP model is a critical example of modified gravity models in cosmology that act as 
alternatives to dark energy. 


The RS and DGP models are 5-dimensional phenomenological models, and so a key issue is 
how to realize such models in 10-dimensional string theory. Some progress has been made. 6- 
dimensional cascading brane-worlds are extensions of the DGP model. 10-dimensional type IIB 
supergravity solutions have been found with the warped geometry that generalizes the RS geometry. 
These models have also been important for building inflationary models in string theory, based on 
the motion of D3 branes in the warped throat [63] [[36] (see the reviews [294] [32] and references 
therein). The action for D3 branes is described by the Dirac-Born-Infeld action and this gives the 
possibility of generating a large non-Gaussianity in the Cosmic Microwave Background temperature 
anisotropies, which can be tested in the future experiments [398] [213] (see the reviews [86] [250)). 

These models reply on the effective 4-dimensional approach to deal with extra dimensions. 
For example, the stabilization mechanism, which is necessary to fix moduli fields in string theory, 
exploits non-perturbative effects and they are often added in the 4-dimensional effective theory. 
Then it is not clear whether the resultant 4-dimensional effective theory is consistent with the 10- 
dimensional equations of motion [108] [109] [238] 251]. Recently there has been a new development 
and it has become possible to calculate all significant contributions to the D3 brane potential in 
the single coherent framework of 10-dimensional supergravity [31] [30] [28] [29]. This will provide us 
with a very interesting bridge between phenomenological brane-world models, where dynamics of 
higher-dimensional gravity is studied in detail, and string theory approaches, where 4D effective 
theory is intensively used. It is crucial to identify the higher-dimensional signature of the models 
in order to test a fundamental theory like string theory. 

In summary, brane-world gravity opens up exciting prospects for subjecting M theory ideas 
to the increasingly stringent tests provided by high-precision astronomical observations. At the 
same time, brane-world models provide a rich arena for probing the geometry and dynamics of the 
gravitational field and its interaction with matter. 
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